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-h ! Abstract 

We estimate analytically the critical coupling separating the weak and the strong coupling 
regime in 4D simplicial quantum gravity to be located at fc§" — 1.3093. By carrying out 
a detailed geometrical analysis of the strong coupling phase we argue that the distribu- 
tion of dynamical triangulations with singular vertices and singular edges, dominating in 
such a regime, is characterized by distinct sub-dominating peaks. The presence of such 
peaks generates volume dependent pseudo-critical points: kf lt (N4 = 32000) ~ 1.25795, 
kf u {N 4 = 48000) ~ 1.26752, k^ u {N 4 = 64000) ~ 1.27466, etc., which appear in good 
agreement with available Monte Carlo data. Under a certain scaling hypothesis we an- 
alytically characterize the (canonical) average value, ci(A4;A:2) =< Nq > /N 4 , and the 
susceptibility, C2(N 4 ; ^2) = (< Nq > — < Nq > 2 )/N 4 , associated with the vertex distribu- 
tion of the 4-D triangulations considered. Again, the resulting analytical expressions are 
found in quite a good agreement with their Monte Carlo counterparts. 
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1 Introduction 



In this article we shall characterize analytically the critical point separating the weak and 
the strong coupling regime in 4_D-simplicial quantum gravity by locating it at k^ lt ~ 
1.3093. The elementary techniques we develop here will allow us to get a rather detailed 
understanding of the geometry and the physics of the strong coupling phase of the theory. 
In particular we will show that the dynamics in such a phase is influenced by the presence 
of peaks in the distribution of singular triangulations. These latter are combinatorial man- 
ifolds characterized by the presence of vertices shared by a number of simplices diverging 
linearly with the volume of the triangulation, and possibly connected by a sub-singular 
edge. The peaks in question are parameterized by the fraction of total volume which is al- 
located around such singular vertices. In order of decreasing entropic relevance, the peaks 
are found, according to a well-defined pattern, at &2 — 1.24465, &2 — 1.2744, A>2 — 1.2938, 
&2 — 1.30746, &2 — 1.31762, k<i ~ 1.32545, etc., asymptotically fading towards the weak 
coupling regime. By exploiting simple entropic arguments drawing from our recent work 
p|, [Q] and by making use of a certain scaling hypothesis, we show how such collection 
of sub-dominating sets of singular triangulations significantly affects the dynamics of the 
transition between weak and strong coupling. We hope that our work offers the possibility 
of making progress in understanding the nature of the transition, one of the major issue 
which controls the validity of dynamical triangulations as the basis of a regularization 
scheme for gravity. Before embarking on this analysis we offer some general motivation 
for such a study. 

1.1 The Model 

Let M be a closed n-dimensional, (n > 2), manifold of given topology. Let Riem(M) and 
Diff(M) respectively denote the space of Riemannian metrics g on M, and the group of 
diffeomorphisms on M. In the continuum formulation of Euclidean quantum gravity one 
attempts to give meaning to a formal path integration over Diff(M) equivalence classes 
of metrics in Riem(M): 



Z(A,G,M)= [ V[g(M)]e- s ^ A ' G '^ (1) 

J Riem(M)/Diff(M) 

where S g [A, G, S] is the Einstein-Hilbert action associated with the Riemannian manifold 
(M,g), viz., 



S g [A,G,E]=A [ d n i^g--^-\ d n ^R (2) 

JM lD7rG Jm 

and T>[g(M)] is some a priori distribution on Riem(M) / Dif f (M) describing the strong 
coupling statistics (A — > 0, G — > oo) of the set of Riemannian manifolds {(M,g)} con- 
sidered. We avoid here discussing well-known specific pathologies in dealing with (|l|) and 
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about which the reader can find abundant literature, and simply recall that in the Dynam- 
ical Triangulations approach to quantum gravity one attempts to give meaning to (|l|) by 
replacing the continuum Riemannian manifold (M, g) with a Piecewise-Linear manifold 
(still denoted by M) endowed with a triangulation T a — * M generated by gluing a (large) 
number of equilateral n-simplices a n . One approximates Riemannian structures by means 
of such triangulated manifolds by using a representative metric where each simplex a n is a 
Euclidean equilateral simplex with sides of length a, (typically we set a = 1). This metric 
is locally Euclidean everywhere on the PL-manifold except near the (n — 2) sub-simplices 
o~ n ~ 2 , (the bones), where the sum of the dihedral angles, 9(a n ), of the incident a n, s is in 
excess (negative curvature) or in defect (positive curvature) with respect to the 2ir flatness 
constraint, the corresponding deficit angle r being defined by r = 2tt — Y^ a n 0(°~ n )- If K n ~ 2 
denotes the (n — 2)-skeleton of T — ► M n , then M n \K n ~ 2 is a flat Riemannian manifold, 
and any point in the interior of an r- simplex a r has a neighborhood homeomorphic to 
B r x C(link(a r )), where B r denotes the ball in W 1 and C(link(a r )) is the cone over the 
link link(a r ), (the product link(a r ) x [0, 1] with link(a r ) x {1} identified to a point). Note 
that for dynamical triangulations the deficit angles are generated by the string of inte- 
gers, the curvature assignments, {q(k)} k =Q 2 1 providing the numbers of top-dimensional 
simplices incident on the N n _ 2 distinct bones, viz., r(i) = 2ir — q(i) arccos(l/n). 

By specializing to this setting the standard Regge calculus, the formal path integration 
([[]) is replaced on a dynamically triangulated PL manifold M, (of fixed topology), by the 
(grand-canonical) partition function j|],(2|],j|] 

Z[k n - 2 ,k n ]= ^e~ k - N - +k ^ N ^ (3) 

where k n -2 and k n are two (running) couplings, the former proportional to the inverse 
gravitational coupling 1/G, while the latter is a linear combination of 1/16ttG and of the 
cosmological constant A. The summation in (||) is extended to the set {T(M)} of all 
distinct dynamical triangulations the PL-manifold M can support, and it is weighted by 
the symmetry factor, Ct, of the triangulation: the order of the automorphism group of 
the graph associated with the triangulation T. Since symmetric triangulations are the 
exception rather then the rule, we shall assume Ct = 1 in the estimates of the partition 
functions below. Thus, in the following we will omit the symmetry factor when writing 
the partition function. 

One can introduce also the canonical partition function defined by 

W(k n . 2 ) eff = e k ^ N -\ (4) 

TeT(N n ) 

where the summation is extended over all distinct dynamical triangulations with given 
N n , (i.e., at fixed volume), of a given PL-manifold M. Finally, we shall consider the 
micro- canonical partition function 

W[N n _ 2 ,b(n,n-2)}= £ 1, (5) 

TeT(N n ;N n - 2 ) 
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where the summation is extended over all distinct dynamical triangulations with given iV, 
and iV n _2, i.e., at fixed volume and fixed average incidence 



of a given PL-manifold M. The micro-canonical partition function is simply the number of 
distinct dynamical triangulations with given volume (oc N n ) and fixed average curvature 
(oc b(n,n — 2)), of a given PL- manifold M. In other words, W[N n _2,b(n,n — 2)] is the 
entropy function for the given set of dynamical triangulations: it provides the discretized 
counterpart of the a priori distribution D[g(M)] describing the strong coupling statistics 
(k n — > 0, k n -2 — ► 0) of the set of Riemannian manifolds {(M, g)}. 

Recently there have been a number of significant advances in 3- and 4-dimensional 
simplicial quantum gravity that fit together in a coherent whole; roughly speaking, these 
results are related to: (i) a deeper understanding of the geometry of n > 3-dimensional 
dynamical triangulations; (ii) the study of simpler models mimicking quite accurately the 
critical structure of simplicial quantum gravity; (iii) more refined computer simulations 
of the phase structure of the 4- dimensional theory. These results imply that both in 
dimension n = 3 and n = 4, simplicial quantum gravity has two geometrically distinct 
phases parameterized by the value of the inverse gravitational coupling /c n _2- In the weak 
coupling phase (large values of 2) we have a dominance of PL-manifolds which collapse 
to branched polymer structures with an Hausdorff dimension du = 2 and an entropy 
exponent (analogous to the string susceptibility of the 2- dimensional theory) 7 = 1/2. 
In this phase the theory has a well defined continuum limit which is independent of any 
fine tuning of the (inverse) gravitational constant k n -2- We are not really interested in 
this continuum limit, even if it exists. The situation is not unusual from the point of 
view of lattice theories. For instance in compact U{\) gauge theories one hits the trivial 
Coulomb phase for all f3 > In the strong coupling phase (small values of A; n _2), 
we have a dominance of crumpled manifolds: the typical configuration sampled by the 
computer simulations is that of a triangulation with a few vertices on which most of 
the top- dimensional simplices are incident, such presence of singular vertices, typically 
connected by a sub-singular edge, seems to be a signature of the strong coupling phase ||. 
Note that the word singular vertex is used in DT theory with a meaning quite different 
from the accepted meaning adopted in PL geometry. What is actually meant is that a 
metric ball around any such a vertex, of radius equal to the given lattice spacing, has a 
volume that grows proportionally to the volume of the whole PL-manifold. This behavior 
indicates that the Hausdorff dimension of the typical triangulation in the strong coupling 
phase is very large if not infinite. There is strong evidence that the transition between 
weak and strong coupling, marked by a critical value kn—2, is of a first order nature in 
the n = 3-dimensional case. In dimension n = 4, the original numerical simulations 
seemed to indicate a second order nature of the transition, a result that invited positive 
speculations on the possibility that simplicial quantum gravity could indeed provide a 
reliable regularization of Euclidean quantum gravity. However, recent and more accurate 
analyses Q of the Monte Carlo simulations seem rather to point toward a first order nature 
of the transition. These results are not definitive since the latent heat at the critical point 
is very small as compared to the 3D-case, so that the question remains open whether 




(6) 
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any further increase in the sophistication of the simulations will definitively establish, 
within the limit of the reached accuracy, the nature of the transition. In any case, all 
recent numerical results jjj strongly indicates that the phase transition in 4-dimensional 
simplicial gravity is associated with the creation of singular geometries. These results have 
put to the fore the basic problems of theory, and in this sense important questions abound: 
What is the geometrical nature of the crumpled phase? What is the mechanism driving 
the transition from the polymer phase to the crumpled phase? Is there any geometrical 
intuition behind the first or higher order nature of the transition? Is it possible to take 
dynamical control over the occurrence of singular vertices and edges which otherwise would 
entropically dominate? 

Some of these questions can be systematically addressed by a detailed but otherwise 
elementary discussion of the geometry of dynamical triangulations along the lines of ||], 
H. This geometrical viewpoint has turned out to be useful and interesting in terms of 
providing an analytical framework within which discuss simplicial quantum gravity while 
at the same time maintaining a strong contact with computer simulations. In this paper 
we discuss in detail how this approach can be further exploited to estimate analytically 
the value of fcf lt , and describing the geometrical properties of the strong coupling phase 
of the 4- dimensional theory. 

2 Large volume asymptotics of the partition functions 

For the convenience of the reader we recall in this section a few basic results of [jl|] that 
we are going to use later on. 

The discretized distribution W[N n _2, b(n, n — 2)] is one of the objects of main interest 
in simplicial quantum gravity, and for n > 3 an exact evaluation of W[N n -2, b(n, n — 2)] 
is an open and very difficult problem. However, since in (||) and (0) one is interested in 
the large volume limit, what really matters, as far as the criticality properties of @ are 
concerned, is the asymptotic behavior of W[N n _2, b(n, n— 2)] for large N n . This makes the 
analysis of W[_/V n _2, b(n, n— 2)] somewhat technically simpler, and according to jjj one can 
actually estimate its leading asymptotics with the relevant sub-leading corrections. If we 
consider an n-dimensional (n > 2) PL-manifold M of given fundamental group TTi(Af)), 
then the distribution W[N n -2, b(n, n — 2)] of distinct dynamical triangulations, with given 
N n ^2 bones and average curvature b = b(n,n — 2), factorizes according to 

W[N n _2,b(n,n-2)]=p c Z v _ 2 (Card{T^} curv ), (7) 

where Pjy r " ^ s ^ e number of possible distinct curvature assignments {g(a)}^Q 2 for tri- 
angulations {To} with iV„_2 bones and given average incidence b(n,n — 2), viz., 

{q(a)}^o 2 + Um?= 2 / M7)fe 2 ^ - , (8) 
while < Card{T^ ] } 

curv > is the average (with respect to the distinct curvature assign- 
ments) of the number of distinct triangulations sharing a common set of curvature as- 
signments, (for details, see section 5.2 pp.102 of |l]]). This factorization allows a rather 
straightforward asymptotic analysis of W[N n ^2, b(n, n — 2)], and in the limit of large N n 
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(theorem 6.4.2, pp. 131 with the expression < Card{Ta^} curv > provided by 



theorem 5.2.1, pp. 106) 

W[N n „ 2 ,b(n,n-2)} 



. e (a„b(n,n-2)+a n -2)N n -2 



'(b-q+l 



il-2n 



_|_ 1)6-9+1 



(b-q) 

-^Kn,^"! ( b(n,n-2) 

n(n + l) " - 



JV»-2 



(9) 



D/2 



N n -2 



r(b)~^ 



The notation here is the following: 

is a topology dependent parameter of no importance for our present purposes (see 
Q for its explicit expression), a n _2 and a n are two constants depending on the dimen- 
sion n, (for instance, for n = 4, a n = — arccos(l/ra)| n= 4, a n _2 = 0); q is the minimum 
incidence order over the bones (typically q = 3); D=dim[Hom(Tri(M),G)] is the topolog- 
ical dimension of the representation variety parameterizing the set of distinct dynamical 
triangulations approximating locally homogeneous G-geometries, (G C SO(n)). Finally, 
m(b) > and r(6) > are two parameters depending on b(n, n — 2) which, together with 
riff > 1, characterize the sub-leading asymptotics of W[N n _2, b(n, n — 2)]. In particular, 
note that 



^\- m (b(n,n~2))Nl /n H] Mn,n - 2) o) D/2 N r(b) 

y n(n + l) " " ' 



(10) 



is the asymptotics associated with (Cor<i{r^} curv ). The remaining part of @ is the 
leading exponential contribution coming from the large N n behavior of the distribution 
Pn™ '2 °^ ^ e P oss ible curvature assignments. This latter term provides the correct be- 
havior of the large volume limit of dynamically triangulated manifolds, an asymptotics 
that matches nicely with the existing Monte Carlo simulations, (e.g., see [jjj, section 7.1, 
pp.160). 



While the exponential asymptotics is basically under control, it must be stressed that 
some of the most delicate aspects of the theory are actually contained in (Card{T^} curv ) . 
Roughly speaking, the set of triangulations {T^} curv can be rather directly interpreted 
as a finite dimensional approximation of the moduli space of constant curvature metrics 
in smooth Riemannian geometry, (the standard example being the parameterization of 
the moduli space of surfaces of genus > 2 with the set of inequivalent constant curvature 
(= —1) metrics admitted by such surfaces; due to rigidity phenomena, such example is to 
be taken with care for n > 3) . With this geometric interpretation in mind we can consider 

\n(Card{T^} curv ) _ _ m(b(n, n - 2))Nn /nH D 
In iV„ ~ lniV n h ~2~ 

as the (formal) covering dimension of such a moduli space. Clearly, whenever m(b) > 
(and tih finite or 0(1nN n )), such covering dimension is singular. This signals the 
fact that in the corresponding range of the parameter b, dynamical triangulations fail to 
approximate in the large volume limit any smooth Riemannian manifold. 



(n) for N n -» oo (11) 
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2.1 The critical incidence 



The parameters m(b), r(b), and njj, characterizing the covering dimension and the large 
volume asymptotics (|To|) are not yet explicitly provided by the analytical results of ||]. 
By exploiting geometrical arguments, one can only prove[|l]] (theorem 5.2.1, pp. 106) 
an existence result to the effect that if n > 3, there is a critical value bo(n), of the 
average incidence b(n, n — 2), to which we can associate a critical value £^"2 °f the- inverse 
gravitational coupling, such that 

m(b) = 0, (12) 

for b(n,n — 2) < &o( n )> whereas 

m(b) > 0, (13) 

for &o("-) < b(n,n — 2). In other words, for b < bo(n) the sub-leading asymptotics in (^) 
is at most polynomial, whereas for b > &o( n ) this asymptotics becomes sub-exponential 
as N n goes to infinity, (note that in the 2-dimensional case (j9|) has always a sub-leading 
polynomial asymptotics). 

This change in the sub- leading asymptotics qualitatively accounts for the jump from 
the strong to the weak coupling phase observed in the real system during Monte Carlo 
simulations. However, the lack of an explicit expression for m{b) hampers a deeper analysis 
of the nature of this transition. In particular, one is interested in the way the parameter 
m(b(n,n — 2)) approaches as b(n, n — 2) — > bo(n), since adequate knowledge in this direc- 
tion would provide the order of the phase transition. It is clear that a first necessary step 
in order to discuss the properties of m(b(n, n — 2)) is to provide a constructive geometrical 
characterization of the critical average incidence bo (n), and not just an existence result. 

As far as the other parameter r(n) is concerned, the situation is on more firm ground. 
r(n) characterizes the sub-leading polynomial asymptotics in the weak coupling phase, 
and recently []|], an analysis of the geometry of dynamical triangulations in this phase has 
provided convincing analytical evidence that r(n) — (2n + 3)/2 + 3 = 1/2. As expected, this 
corresponds to a dominance, in the weak coupling phase, of branched polymers structures. 



2.2 Canonical Averages and the Curvature Susceptibility 



If we consider the weighted distribution W[N n _2, b(n, n—2)] exp[£; n _2-/V„_2]! characterizing 
the canonical partition function ([|), then it is straightforward to check that, as N n — > 
00, this distribution is strongly peaked around triangulations with an average incidence 
b*(n,n — 2;k n _2) given by (see Q, eqn. (6.39), pp. 134) 



b*(n, n — 2; k 



n-2; 



3( 



A(k 



n-2) 



A(k n _ 



where for notational convenience we have set 
A(k n - 2 ) = 



27 k „ . ,27 



Y e + 1 + V ( y e fc ™-2 + l) 2 



27 



,kn-2 



+ 1 



27 

(_ e fcn-2 + 1)2 



1/3 



1/3 



1. 



(14) 



(15) 



7 



This remark allows to compute, via a uniform Laplace estimation, the large volume asymp- 
totics of the canonical partition function W(h n -%) e ff = J2reT(N n ) e kn ~ 2Nn - 2 . A discussion 
of this asymptotics at the various orders is rather delicate and the reader can find the de- 
tails in [0, (chapter 6, theorem 6.6.1). For our purposes it is sufficient to consider the 
leading order expression which can be readily obtained, starting from the micro-canonical 
partition function, by means of a standard saddle point evaluation, viz. 

W(N n , k n - 2 ) eff = c n ( {A [ k ;- 2) \ 2) ) n N^+ D ^- ^l-rnmn^k^N^] . 



where c n is a scaling factor not depending from h 2 . 



As the inverse gravitational coupling h n - 2 varies, the average curvature correspond- 
ent 



ingly changes according to (H). It follows that there is a well defined critical value, k CTtS 



solution of the equation 



W = 3 ^rf^ a) i )' (17) 

A(k n - 2 ) ~ 1 



for which b*(n,n — 2; h n - 2 ) = £>o( n )> where &o( n ) is the critical average incidence, (see 
(|TJ) and ( p^) ). This fc£"* 2 describes the transition between the strong coupling phase 
(k-n-2 < ^n-l) an d the weak coupling phase {h n _ 2 > associated with the two 

distinct sub- leading asymptotics regimes of fllBp. 



The explicit geometric characterization of &o ( n ) an d the evaluation of the corresponding 
critical value in the 4-dimensional case are among the most important issues that we 
discuss in this paper. In order to compare the geometrical results we obtain with the data 
coming from recent Monte Carlo simulation it will be useful to have handy the expressions 
of the free energy In W e ff(N±, hi), of the (canonical) average of the number of bones 
< N 2 >= dlnW e ff(N4,h 2 )/dh 2 , and of the associated curvature- curvature correlator 
[(iVf) - (N 2 ) 2 } = d 2 lnW eff (N 4 ,k 2 )/dh 2 2 . 

The large volume asymptotics of the (canonical) free energy is readily obtained from 
(|l6|); by setting n = 4 and discarding the inessential constant terms we get (in the saddle- 
point approximation used above) 

In W(N, h 2 ) eff = IOA4 In A ^ + 2 _ mih^N 1 /^ . (18) 



The canonical average < N 2 > follows from differentiating ([l^) with respect to the 
inverse gravitational coupling h 2 , 

/AM - in Al( - k ^ TV A A/n H dm{h 2 ) 

{N2) ~ W A(h 2 ) + 2 Ni ~ N * (19) 
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where we have set 

dA(k 2 ) 



A 



dko 



9 

-e 
2 



1 + 



f e fe 2 + 1 



fe fc 2 + l)2_i 



&2 



e K2 + 1 



y e* 2 + 1 + 



fe^ + 1 



27 

— e fc 2 + 1) 2 -1 
' 2 7 



-2/3 



+ 



27 

(fiefca +1)2-1 

v 2 ; 



-2/3 



(20) 



Note that for a 4-dimensional PL manifold M of Euler characteristic %(M), we have 
A^o = 2 — X- Thus, from (JV2) we immediately get the expression for the first 

normalized cumulant of the distribution of the number of vertices of the triangulation, 
viz. 

r i/n H -idm(k 2 ) 



(Ar , , . (No) 5A 1 (k 2 ) 

C 1 (iV4^2) = ^- = IWT ^-l-iV 4 



dk 2 



(21) 



This is a typical quantity monitored in Monte Carlo simulation, and later on we will 



discuss how (21) actually compares with respect existing numerical data. 
Finally, the curvature susceptibility 



(Ni) - (N 2 



N 4 



J_ <P lnW eff (N 4 ,k 2 ) 
N 4 



dk 2 2 



(22) 



is explicitly computed as 
4c 2 (iY 4 ;/c2) = 



(Nj) - (jV 2 ) 2 

(A(k 2 ) + 2)A 2 (k 2 )-A 1 (k 2 ) 2 _ y nH -i d 2 m(k 2 



{A(k 2 ) + 2Y dk 2 2 ' (23) 

where c 2 (N 4 ;k 2 ) = ^ol_jj^o) - g se cond normalized cumulant of the distribution of 
the number of vertices of the triangulation, and where we have set 



d 2 A(k 2 ) 



A 2 {k 2 ) = n / 2 = A l {k 2 ) + 
2 

243 r 2k 9 . 



+ 



dk 2 - 
-(- 



[(^ e k2 + 1)2 _ l]3/2 

243 r 2k 9 . 

4_; 

[(2* e k2 + 1)2 _ 1]3 /2 



T« k2 + 1 



(ye^ + 1) 2 -1 



^ + 1 + ^(^ + 1)2-1 



-2/3 



-2/3 



81 



2k 2 



f ^ + 1 + ^(^ + 1)2-1 



81 



2fc 2 



(y e fc 2 + 1)2-1 



-5/3 



-5/3 



1 + 



f e fc2 + 1 



(fe*> + 1)2-1 
f e fc2 + 1 

(27 e fe 2 + 1)2 _ ! 



(24) 
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2.3 A scaling hypothesis 

Clearly the above expressions for c\(N±; k 2 ) and c 2 (N^; k 2 ) are useless if we do not specify 
how m(k 2 ) depends on the inverse gravitational coupling k 2 . Since according to theorem 
5.2.1 of 0], m{k 2 ) — » as b(k 2 ) approaches a critical incidence 6o(4), (henceforth denoted 
by 6o); the simplest hypothesis we can make is that, for (b(k 2 ) — 60) — * + , m(k 2 ) scales 
to zero according to a power law given by 

v b(k 2 ) b 

where < v < 1 is a critical exponent to be determined; (the factor 1/u is inserted for 
later convenience). 

The expression ( |l6| ) for the canonical partition function W[N n _2, b(n, n — 2)] is a large 
volume (N4) asymptotics evaluated at fixed volume. It contains a non- trivial subleading 
asymptotics governed by m(k 2 ). The net effect of this subleading term is clearly visible 
in the expressions (|2l| ) and ( |23|) of the two cumulants, and shows that in order to capture 
the behavior of ci(N4;k 2 ) and c 2 (N4;k 2 ) as k 2 — ► k 2 cnt , ( |25| ) is not sufficient. It must be 
combined with a finite scaling hypothesis telling us how m(&2) scales with the volume N4, 
as b(k 2 ) —>■ bo- From the asymptotics ([l6|), and the expression ( ^T|) for the first cumulant, 
it easily follows that the simplest, if not the most natural, hypothesis we can make is to 
assume that m(k 2 ) scales asymptotically with the volume according to 

hm I -L- --V-i.jVV^-^1, (26) 
Af 4 ^oo, o{k 2 ) b 

where according to theorem 5.2.1 of [|lj uh > 1- This implies that, when b(k 2 ) — » 60, 
m{k 2 ) scales as 

i/(n g -l) 

m(fe) ~ iV 4 . (27) 

Together with m{k 2 )N^ nH — ► as 6(^2) — ► &o> (p6|) yields (1/njf) < ^ < 1, a finite size 
scaling relation connecting the critical exponents u and n#. 
Introducing this ansatz in (^Tj) and ( |23"|) we explicitly obtain 

, , . (iVo) 5^(fc 2 ) 1 ,„ a . 



and 



C2W;fe2) _ MkJ^l! (29) 



iV 4 

5 (A(fc 2 ) + 2)A 2 {k 2 ) - A^ktf 1 A 2 (k 2 )A(k 2 ) - 2A 1 (k 2 ) 
2 (^(fc 2 ) + 2) 2 12 A*(k 2 ) 

>-l| A\{k 2 ) A{k 2 )-l 1 -' 
36 A 4 (fc 2 ) 3A(& 2 ) bo 



+ 
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Note that in this latter expression the only undetermined parameters are the critical 
exponent v and the critical incidence bo . In the following paragraphs we provide an explicit 
geometric characterization of such 6q- The only remaining unknown quantity is then v. 
Anticipating the conclusion of the paper, it turns out that it is possible the choose a value 
of v (~ 0.94) which leads to quite a good agreement between (p^)-(|29|) and the available 
Monte Carlo data for the distributions of these two cumulants. 

3 The Geometry of the strong coupling phase 

A PL manifold endowed with a dynamical triangulation is a particular example of an 
Alexandrov space, i.e. finite dimensional, inner metric space with a lower curvature 
bound in distance-comparison sense, (a brief introduction with the relevant references 
can be found in Q, section 3.2). The natural topology specifying in which sense dynami- 
cal triangulations approximate Riemannian manifolds is associated with an Hausdorff-like 
distance introduced by Gromov ||, and which is a direct generalization of the classical 
Hausdorff distance between compact subsets of a metric space. The role of this topology 
stems from the fact there are many geometric constructions in dynamical triangulations 
theory that are close in Gromov-Hausdorff topology, but not in smooth Riemannian ge- 
ometry. In Q we proved that every Riemannian manifold (of bounded geometry) can be 
uniformly approximated in this topology by dynamical triangulations, ([|]], section 3.3, 
Th.3.3.1); the converse result, namely if every dynamical triangulation approximates, as 
the number of simplices goes to oo, a n-dimensional Riemannian manifold, is deeply tied 
to understanding the structure of the thermo-dynamical behavior of the large volume limit 
of the set of possible dynamical triangulations. It is interesting to note that the geometry 
of the set of all possible dynamical triangulations of a manifold of given topology is a sub- 
ject of which little is actually known, even in dimension two. Recently, in a remarkable 
paper [jlC| 1, W. Thurston has shed some light in the two- dimensional case by showing that 
the space of triangulations (of positive curvature) has the rich geometric structure of a 
complex hyperbolic manifold. We do not need to reach, in this work, such a level of so- 
phistication and in order to determine the critical average incidence, bo, we discuss mostly 
the kinematical properties of the space of all possible dynamical triangulations admitted 
by an n- dimensional PL manifold M of given topology. 

Let (M, T a ) be a dynamically triangulated manifold, then the /- vector of the trian- 
gulation is the string of integers (N (T a ), Ni(T a ), . . . ,N n (T a )), where Ni(T a ) G N is the 
number of i-dimensional sub- simplices a % of T a . This vector is constrained by the Dehn- 
Sommerville relations 

n 

5>im(T) = X (T), (30) 

i=0 

i'-' v*^-^ ^ <3i> 



ii 



if n is even, and 1 < k < n/2. Whereas if n is odd 

yC-lY-, ,, Ni(T) = 0, (32) 

^ ' (i-2k + iy.2k\ K ' y 1 

i=2k v ' 

with 1 < k < (n — l)/2, and where x{T) is the Euler- Poincare characteristic of T. It 
is easily verified that the relations (|30|) , ([H]), (|32|) leave — 1, (n even) or i(n — 1) 
unknown quantities among the n ratios Ni/Nq, . . . , N h /Nq, [fllf] . Thus, in dimension 
n = 2, 3, 4, the datum of N n , and of the number of bones iV n _2, fixes, through the Dehn- 
Sommerville relations all the remaining iVj(T). These extremely simple and perhaps even 
naive-sounding remarks turn out to be quite powerful in providing information on the 
global metrical properties of the underlying PL-manifold. Not only, as is obvious, on the 
volume (oc N n (T a )), and on the average curvature (oc ^n(n + l)N n (T a )/N n _2(T a )), but, 
corroborated by a few more elementary facts, also on the genesis of singular vertices and 
edges. 

An elementary but geometrically significant result of this type is provided by the range 
of variation of the possible average incidence 6(n, n — 2). One gets (see|l|] (lemma 2.1.1)) 

Lemma 1 Let T a — > M n a triangulation of a closed n- dimensional PL- manifold M , with 
2 < n < 4 ; then for N n (T a ) — > oo, we get 

(i) For n = 2: 

6(2,0) =6; (33) 

(ii) for n = 3: 

~ < 6(3, 1) < 6; (34) 

(Hi) for n = 4; 

4 < 6(4, 2) < 5. (35) 



The 2-dimensional case as well as the upper bounds for n = 3 and n = 4 are well-known 
trivial consequences of the Dehn- Sommerville relations. The lower bounds 6(3,1) > 9/2 
and 6(4, 2) > 4 are instead related to a rather sophisticated set of results proved by 
Walkup[12] in the sixties concerning the proof of some conjectures for 3- and 4-dimensional 
PL manifolds; (apparently, these results went unnoticed by researchers in simplicial quan- 
tum gravity). Walkup's theorems have important implications for understanding the geom- 
etry both of the strong and of the weak coupling phase of simplicial gravity. In dimension 
n = 3, we have 111 

Theorem 1 There exists a triangulation T — > S 3 of the 3-sphere § 3 with Nq vertices and 
Ni edges if and only if Nq > 5 and 

4jV -10> Nl > N ^~ 1 \ (36) 

Moreover T is a triangulation of § 3 satisfying Aq = 4Nq — 10 if and only if T is a stacked 
sphere, whereas T is a triangulation of S 3 satisfying Aq = jV °(^°~ 1 ) if an d on \y if T is a 
2-neighborly triangulation, namely if every pair of vertices is connected by an edge. 
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A stacked sphere (S™, T) is a triangulation T —* § n of a sphere which can be constructed 
from the boundary da n+1 ~ § n of a simplex <r n+1 by successive adding of pyramids over 
some facets. More explicitly, the boundary complex of any abstract (n + l)-simplex a n+1 
is by definition a stacked sphere, and if T is a stacked sphere and a n is any n-simplex of T, 
then T is a stacked sphere if T is any complex obtained by T by removing <r n and adding 
the join of the boundary da n with a new vertex distinct from the vertices of T. Note also 
that a triangulated PL-manifold is called /c-neighborly if 

= WK^y. (37> 

We are referring explicitly to 3- and 4-spheres S n , because the majority of Monte Carlo 
simulations have been carried out in these cases (for a recent discussion of more general 
topologies, see[|n|). However, it must be stressed that the above definitions, as well as 
Walkup's theorems, can be naturally extended (with suitable modifications |l2|]) to any 
n-dimensional PL manifolds M. Note in particular that every triangulable 3- manifold M 
can be triangulated so that the closed star of some edge contains all the vertices and every 
pair of vertices is connected by an edge. 

In dimension n = 4 we have a somewhat weaker characterization of the possible set of 
triangulations: 

Theorem 2 IfT^Misa triangulation of a closed connected A-manifold, then 

Ni(T) > 5N (T) - y X (T), (38) 
and equality holds if and only if (M, T) is a stacked sphere. 

Note that actually one has a stronger statement in the sense that equality in ( |3"8| ) holds if 
and only is all vertex links in the 4-manifold M are stacked 3-spheres. 

Contrary to what happens for 3-manifolds, 2-neighborly triangulations (i.e., trian- 
gulations where every pair of vertices is connected by an edge), are not generic for 4- 



dimensional PL manifolds, and as matter of fact, the above theorem immediately implies [14] 
that for any such (M, T) 

N (T)(N (T) - 11) > -15 X (M), (39) 

where the equality implies that (M, T) is 2-neighborly. Thus, the equality is not possible 
for large and arbitrary values of Nq(T), but, (depending on topology) [[U|] , only in the cases 
N (T) = 0, N = 5, iV = 6, or iV = 11 mod 15. 

Even if 2-neighborly triangulations are not generic, one can easily construct voluminous 
(i.e., with A^(T) arbitrarily large) triangulations of the 4-sphere where all vertices but two 
are connected by an edge. In order to realize such triangulations, consider a 2-neighborly 
triangulation T(3) of the 3-sphere S 3 with /-vector [JV (T(3)), Aq(T(3)), iV 2 (T(3)), JV 3 (T(3))]. 
If we take the Cone, C(§ 3 ), on such (§ 3 ,T(3)), viz., the product § 3 x [0,a] with S 3 x {a} 
identified to a point, then we get a triangulation of a 4-dimensional ball B 4 with /-vector 
given by 

f(B A ) = (iVo(T(3)) + l,Af 1 (T(3)) + iVo(T(3)), 

iV 2 (T(3)) + JVi(r(3)), N 3 (T(3)) + N 2 (T(3)), N 3 (T(3))). (40) 
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By gluing two copies of such a cone C(§ 3 ) along their isometric boundary <9C(8 3 ) ~ S 3 , 
we get a triangulation of the 4-sphere § 4 with /-vector 

f(S 4 ) = (JV (T(3)) + 2,JVi(r(3)) + 2JVo(T(3)), 

N 2 (T(3)) + 2N 1 (T(3)),N 3 (T(3)) + 27V 2 (T(3)), 2iV 3 (T(3))). (41) 

It is trivially checked that corresponding to such a triangulation we get 

34W i\r 



(42) 



where the —1 accounts for the missing edge between the two cone- vertices in C(§ 3 ) U53 
C(S 3 ). 

When applied to simplicial quantum gravity, the existence of such 2- neighborly (or 
almost 2-neighborly) triangulations implies that there are dynamical triangulations of 
the n-sphere 8™, n = 3, n = 4, where all vertices are singular. Corresponding to such 
configurations we have that b(n,n — 2)| n= 3 = 6, and b(n,n — 2)| n= 4 = 5. Thus, not 
surprisingly, for such triangulations the kinematical upper bound for the average incidence 
b(n, n — 2) is attained. However, it is important to stress that such extremely singular 
configurations do not saturate the set of possible configurations for which b max (n,n- 2) 
is reached. From the Dehn-Sommerville relations one immediately gets 

6(n,n-2)| n . 3 = 6.^-, (43) 

and 

6(n,ra-2)L =4 = 10- — -4 — , (44) 

v ' n 2N A + 2N - 2x(T) v ; 

which, together with the obvious relation Ni < Nq(Nq — l)/2, implies that in order to 
attain the upper kinematical bounds b max (n, n — 2)| n= 3 = 6 and b max (n 1 n — 2)| n= 4 = 5 it 
is sufficient that 



N (T) = 0[N n (T) a }, (45) 

with 1/2 < a < 1. Note that 2-neighborly or almost 2- neighborly triangulations corre- 
spond to a = 1/2. 



3.1 Singular Stacked Spheres 

It should be stressed that the presence of singular vertices can occur also for b(n, n — 2) = 
b m in(n, n — 2), i.e., for stacked spheres. In other words, singular vertices are not kinemat- 
ically forbidden by the geometry of the triangulations. Their suppression or enhancement 
in the different phases of simplicial quantum gravity is rather related to the relative abun- 
dance, with respect to the totality of possible triangulations, of the number of distinct 
triangulations with singular vertices as b(n, n — 2) varies. In other words, it is an en- 
tropic phenomenon as clearly suggested by S. Catterall, G. Thorleifsson, J. Kogut, and R. 
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Renken For definiteness, we can describe a concrete construction of a singular stacked 
sphere. It amounts to gluing a 4-dimensional ball B 4 bounded by a stacked 3-sphere S 3 
with a cone over such an S 3 . 

Consider a 3-dimensional stacked sphere § 3 . According to one of Walkup's theorems, 
such an § 3 is the boundary of a 4-dimensional ball B 4 with a tree-like structure and 
corresponds to a triangulation with /-vector 

f(B 4 ) = (N (S 3 ), A^S 3 ), N 2 (S 3 ), N 3 (S 3 ) + N 3 (B 4 ), N 4 (B 4 )), (46) 

where Ni(S 3 ), z = 0,1,2,3 is the /-vector of the boundary stacked sphere and N 3 (B 4 ) is 
the number of a 3 in the interior, B 4 , of B 4 . Note that if we take the cone, C(5 3 ), over the 
boundary stacked sphere, we get another triangulation of the 4-dimensional ball, B 4 ing , 
whose boundary is again isometric to the given S 3 , but whose interior contains a (unique) 
singular vertex. The /-vector of such a triangulation is 

f(B% ng ) = (l+N (S 3 ),N 1 (S 3 )+N (S 3 ),N 2 (S 3 )+N 1 (S 3 ),N 3 (S 3 )+N 2 (S 3 ),^ (47) 

Gluing these two triangulated balls B to B 4 ing through their common boundary § 3 , we 
get a triangulation of the 4-sphere § 4 ~ B 4 \J s s C(dB 4 ), with /-vector 

N = N (S 3 ) + 1 

Ni = N^S 3 ) + N (S 3 ) 

N 2 = N 2 {S 3 ) + iVi(5 3 ) 

N 3 = N 3 (S 3 ) + N 3 (B 4 ) + N 2 (S 3 ) 

iV 4 = N 4 {B 4 ) +N 3 (S 3 ). (48) 

Since S 3 is a stacked sphere, we have 4N 3 (S 3 ) = 3N\(S 3 ) which, together with N 2 (S 3 ) = 
2N 3 {S 3 ), implies 

N 2 = l ^N 3 {S 3 ). (49) 

From 2N 3 = hN± and the Euler relation for the 4-dimensional ball B 4 , (with x(B 4 ) = 1), 
we immediately get N^{B 4 ) = ^N 3 (S 3 ) — |, which implies 

iV 4 = ^V 3 0S 3 )-! (50) 

Thus, for N 3 (S 3 ) — > oo we get a voluminous triangulation (N^ — ► oo) of § 4 with average 
incidence 

N = -N ( S 3 ) — - 
6(n,n-2)L_ 4 = 10- — 3 in 31 / ON 3 ->jv,^oo 4, (51) 

which shows that § 4 ~ B 4 U S 3 C{dB 4 ) is a stacked sphere with a singular vertex (the apex 
of the cone C(dB 4 )). An even simpler construction suffices to prove an analogous result 
in the 3-dimensional case. 



15 



As mentioned in the introductory remarks, stacked spheres are relevant in provid- 
ing the geometrical rationale for the prevalence of branched polymer structures in the 
weak coupling phase of simplicial quantum gravity. As a matter of fact ||],it is their 
tree-like structure that accounts for the kinematical possibility of polymerization. How- 
ever, the existence of stacked spheres with singular vertices, shows that the dynamical 
onset of polymerization is not just a consequence of the geometry of triangulated man- 
ifolds as b(n,n — 2) — ► b m i n (n,n — 2). On the kinematical side we may have, in the 
configuration space, extremal cases such as the 2-neighborly triangulations occurring for 
b(n,n — 2) — > b max (n,n—2) or the singular stacked spheres for b(n, n — 2) — > b m i n (n,n — 2). 
Monte Carlo simulations do confirm that such configurations are not generic. Near 
b max (n,n — 2) we generically sample singular triangulations with just a few singular 
vertices ]|]. Similarly, as b(n, n — 2) — ► b m i n (n, n — 2) the dominant configurations sampled 
correspond to stacked spheres without singular vertices. The mechanism for understand- 
ing the dynamical prevalence of such configurations over the other configurations which 
are kinematically possible is simply related to the fact that, with respect to the counting 
measure, distinct dynamical triangulations are not equally probable as a function of the 
average incidence b(n,n — 2). In order to discuss this point we need to exploit a few 
elementary facts related to the geometry of the ergodic moves used in simplicial quantum 
gravity. 



3.2 Ergodic moves and the onset of criticality 

The (k,l) moves [jlq] in 3 and 4 dimension are a well known set of elementary surgery 
operations (related to the Pachner moves [|l6| ) which allow to construct all triangulations 
of a PL- manifolds starting from a given triangulation. Roughly speaking, the generic 
(k, I) move consists in cutting out a sub-complex made up of fc-dimensional simplices a k 
and replacing it with a complex of I- dimensional simplices a 1 with the same boundary. 
Note that k + I = n + 2. We are interested in discussing how a finite set of such moves 
generate the /-vector of voluminous triangulations of the n- sphere § n , (n = 3, 4), starting 
from the standard /-vector of the simplex da n+1 ~ S n . For n = 3, the relevant moves 
are the (1,4) move (barycentric subdivision), the (2,3) move (triangle to link exchange) 
and their inverses. For n = 4, since the flip move (3, 3) does not alter the distribution of 
the number iVj(T) of simplices, the /-vector of the sphere is generated by the moves (1,5) 
(barycentric subdivision) and (2,4) (two-four exchange) and their inverses. Following^], 
(with a slight change in notation), we denote by Pk,i{n) the number of moves of type 
(k,l) in dimension n, and introduce the balance variables (n = 3): x\ = (-Pi,4 — P4,i), 
%2 = (-P2.3 - ^3,2); and (n = 4): y x = (P 1)5 - P 5)1 ), y 2 = (P 2) 4 ~ ^4,2) ■ In terms of such 
quantities we can easily characterize the string of integers {-/Vj}, i = 0, . . . ,n, which are 
possible /-vectors of triangulated § n . 

For n = 3, we get for /(§ 3 ) = (N (S 3 ), N^S 3 ), N 2 (S 3 ), N 3 (S 3 )): 

N (S 3 ) = 5 + xi, 

N^S 3 ) = 10 + 4xi+X2, 

N 2 (S 3 ) = 10 + 6xi + 2x 2 , 

N 3 (S 3 ) = 5 + 3xi+x 2 , (52) 



16 



whereas, for n = 4 we have for /(§ 4 ) = (N (S 4 ), iVi(<S 4 ), N 2 (S 4 ), N 3 (S 4 ), N A (S 4 )): 

N (S 4 ) = 6 + y h 

N^S 4 ) = l5 + 5 yi +y 2 , 

N 2 (S 4 ) = 20 + 10yi + % 2 , 

N 3 (S 4 ) = 15 + 10yi + 5y 2 , 

N 4 (S 4 ) = 6 + 4 Vl + 2y 2 . (53) 

Note that not all /(§ n ) obtained in this way are actual /- vectors of triangulated § n . This 
is a consequence of the fact that the above relations between the {iVj} and the variables 
Pk,l( n ) are equivalent to the Dehn-Sommerville constraints. And these latter are known 
to be necessary but not sufficient conditions in characterizing the possible /-vectors of a 
triangulated manifold; (sufficient conditions have been conjectured by R. Stanley [pT7f|-see[]l|1 
for a brief discussion of this point). 

Walkup's theorems imply the following kinematical bounds on the variables Xi, yi, 
(i = l,2): 

xi > 0, 

VI > 0, (54) 
(both from the obvious condition No(S n ) > n + 2); 

x 2 > 0, 

V2 > 0, (55) 



(the former from N^S 3 ) > 4N (S 3 ) - 10; the latter from N^S 4 ) > 5N (S 4 ) - ^- X (S 
with X (S 4 ) = 2); 

x\ + x\ — 2x 2 > 0, 

y 2 i+yi-2y 2 > 0, (56) 

(both from Ni(S n ) < No(S n )(No(S n ) — l)/2). Finally, one can express the average inci- 
dence b(n, n — 2) as a function of Xi and yi, so as to obtain 

b(n,n-2)| n=3 = 6- 5 + 3 ; 1 + * 2 , (57) 
10 + 4xi + x 2 



and 



b(n,n-2 n=4 = 10- — — — . (58) 
20 + Wyi + 4y 2 



It is also interesting to discuss in terms of the variables X{ and yi , the average incidence 
of the top-dimensional simplices a n on the vertices a of the triangulations considered. A 
straightforward computation provides 

3W = ^E^°) = (n+l)^, (59) 
No {a0} N 
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yielding 



Q{n)\ n=3 = 4- 



5 + 3xi + x 2 



(60) 



5 + xi 



and 



Q(n)\ n=4 = 5 • 



6 + 4j/i + 2y 2 
6 + yi 



(61) 



As expected, Q(n) is not bounded above: when the move (2,3) (for n = 3), or (2,4) 
(for n = 4) dominates, i.e., near the b max (n,n — 2) kinematical boundary, Q(n) — > oo 
as iV n — > oo. One may wonder if this unboundedness is related to the unboundedness of 
the Einstein-Hilbert action, the answer is most likely no. It is certainly reasonable to put 
restrictions on Q(n) in the search of a continuum limit of the theory, and this may change 
the phase structure of the theory. But if it does, it is just an illustration of the fact that 
this particular part of the phase diagram has no relevance for a genuine continuum limit. 
There should be a reasonable universality. This is nicely illustrated in 2-D dynamical 
triangulation theory where any restriction (except the strict flatness constraint q(i) = 6) 
leads to 2D-gravity. 

The above elementary remarks are a trivial restatement of the well known fact that the 
moves (1,4) and (1,5) (the barycentric subdivision) drive the system into the elongated 
phase, whereas the moves (2,3) and (2,4) drive to the crumpled phase. The crumpling 
transition occur as soon as singular vertices are statistically enhanced by the presence of 
enough (2,4) moves with respect to (1,5), (for n = 3 this enhancement is generated by 
the dominance of (2,3) moves with respect to (1,4) moves). 

3.3 The genesis of singular vertices: S 4 V 

In order to characterize the onset of crumpling we describe the /- vector of the generic 
triangulation of S" in a way that clearly shows the mechanism of formation of singular 
vertices. Such a description is obtained by gluing a triangulated ball B n to the cone over 
its boundary dB n ~ § n_1 . Thus, by referring to the 4- dimensional case for definiteness, 
we consider B> 4 V ~ B 4 U53 C(dB 4 ), (sv for singular vertex). Note that any triangulation 
of S 4 can be factorized in this way (since C(dB 4 ) and dB 4 are the star and the link of a 
vertex, respectively), and we have 



namely when the cone C(dB ) contains a number of top- dimensional simplices growing 
linearly with the volume of the whole manifold. 

It is easily checked that the /-vector of S 4 S ~ B A U53 C{dB A ) is given by 



N 4 = N 4 {B 4 ) + N A (C(dB 4 )). 



(62) 



The triangulation is singular as soon as we have 



N A {B 4 ) oc N A {C{dB 4 )), 



(63) 



iVo = iV (<S 3 ) + l + iVoCB 4 ) 
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iVi = ^(S 3 ) + N {S 3 ) + N^B 4 ) 

N 2 = N 2 (S 3 ) + N^S 3 ) + N 2 (B 4 ) 

N 3 = N 3 (S 3 ) + N 2 (S 3 ) + N 3 (B 4 ) 

N 4 = N 3 (S 3 ) + N 4 (B 4 ), (64) 

where Ni(S 3 ) denotes the /-vector of the boundary d(B 4 ) ~ § 3 of the triangulated ball 
B 4 , and Ni{B 4 ) is the /-vector of the interior of B 4 . The Dehn-Sommerville relations for 
§ 4 and for § 3 constrain Ni(B 4 ) and N k (S 3 ) according to 

N (B 4 ) - Ni(B 4 ) + N 2 (B 4 ) - N 3 (B 4 ) + N 4 {B 4 ) = 1 

2N l {B 4 ) - 3N 2 {B 4 ) + AN 3 {B 4 ) - 5N 4 (B 4 ) + N (S 3 ) = 

2N 3 (B 4 ) + N 3 (S 3 ) = hN A {B 4 ). (65) 

The average incidence 6(4, 2) of such triangulated § 4 can be easily computed in terms 
of the /-vectors Ni(B 4 ) and Nk(S 3 ) according to 

6f42 ) = io 4b (3, l)N 3 {S 3 ) + 26(3, l)[iVo(^ 4 ) - ^Vi(^ 4 ) + N 2 (B 4 )} - 26(3, 1) 

[66(3, 1) + 18]iV 3 (5 3 ) + 36(3, l)Af 2 (5 4 ) 

where 6(3,1) = 6[A r 3(5 3 )/A/'i(5 3 )] is the average incidence of dB 4 ~ § 3 . The presence of 
a singular vertex corresponds to 

MB 4 ) nm (( , 7 , 

and it is easily verified that under such condition 6(4, 2) is an increasing function of 6(3, 1). 
This remark implies that singular triangulations with the smallest possible 6(4, 2) are to be 
found corresponding to 6(3, 1) = 6(3, \) m in = 9/2. 

We have already seen an example of such a triangulation in the previous section, one for 
which the lowest kinematically possible incidence, 6(2,4) = 4, is attained. However such 
examples are not generic. They correspond to assuming y 2 = 0, (or more generally, they 
still occur if one interprets the right hand side as y 2 = O(l)), and the singular vertex is 
not stable under (1,5) moves. Eventually by performing enough barycentric subdivisions 
the initial singular vertex is smoothed out. Explicitly, assume that we start our chain of 
barycentric subdivisions on an a S> 4 V ~ B 4 Ug3 C(dB 4 ) with a given value of N4, say A^(0). 
Denote by § 4 „(0) this initial triangulation. Note that at this initial step 

N,(B\0)) = ±N 4 (C(dB 4 (0))), (68) 

(see the previous section). If we carry out a (1,5) move on each 4-simplex of S 4 ^(0), we 
get a triangulation of § 4 still of the form S 4 V ~ B 4 U53 C(dB 4 ), which we denote S 4 „(l). 
For such triangulation we have 

N 4 (C(dB 4 (l))) = 4 • N 4 (C(dB 4 (0))), (69) 
N 4 (B 4 (1)) = 5 • N 4 (B 4 (0)) + N 4 (C(dB 4 (0))). (70) 
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Now proceed by induction, noticing that if at each step we carry out a barycentric subdi- 
vision of each 4-simplex of the S 4 ^ generated at the previous step, we still get a 4-sphere, 
Sg V (k), triangulated according to S 4 „ ~ B 4 U s a C(dB 4 ), and such that, at the k-th step, 

N 4 (C(dB 4 (k))) = 4 k ■ N 4 (C(dB 4 (0))), (71) 



k-l 

N 4 {B 4 {k)) = 5 k ■ iV 4 (£ 4 (0)) + N 4 (C(dB 4 (0))) ]T 5*" 1 • 4*" 1 . (72) 

i=i 

Thus, as k grows (corresponding to y\ — > +oo, y 2 = 0(1)), N 4 (B 4 (k)) largely dominates 
over N A {C{dB 4 {k))): 

N A (B 4 (k)) » yi ^ +OD N 4 (C(dB 4 (k))), (73) 

and the resulting triangulation of S 4 is no longer singular. In this sense, the dominance of 
the (1,5) move naturally yields regular stacked sphere and thus for branched polymers. 



Discarding these particular examples of unstable triangulated spheres with a singular 
vertex, we can easily characterize the smallest 6(4, 2) corresponding to generic singular 
triangulations, namely triangulations generated in the large volume limit as y\ — > oo and 
y 2 — ► oo, and whose singular vertices are stable under the action of the (k, I) moves (at a 
fixed ratio j^-). Let us start by noticing that corresponding to 6(3, 1) = 9/2, the expression 
(|6G|) for the average incidence reduces to 

b(4 2) = 10 6Nl{S3) + A[N ° { ^ ] - Nl{ ^ ] + N2{ ^ )] (74) 
1 ' ' 15JVi(5 3 ) + 6N 2 {B 4 ) ' 1 ' 

where, in the numerator, we have discarded terms which are o(l), thus irrelevant in the 
large volume limit. 

Since S 3 is stacked, the integers N 3 (S 3 ) and N 4 (S 3 ) are related by 4N 3 (S 3 ) = SN^S 3 )- 
10, which implies that SN^S 3 ) = 10(mod 4). Thus, N^S 3 ) must be an integer multiple 
of 4 up to an error term which goes to zero, with increasing N 4 (S 3 ), as 10/N 4 (S 3 ). More 
explicitly, and referring to the expression of the /-vector of § 3 in terms of the balance 



variables x\ G N and x 2 G N introduced in section |3.2| , we get the following components: 

N (S 3 ) = 5 + xi, 
iVi(S 3 ) = 10 + 4xi, 
N 2 (S 3 ) = 10 + 6xi, 

N 3 (S 3 ) = 5 + 3xi, (75) 
since corresponding to a stacked S 3 we have x 2 = 0, (see (|52|)). 



The congruence properties just established for the /-vector of a stacked 3-sphere sug- 
gest to parameterize both N 2 (B 4 ) and Nq(B 4 ) — Ni(B 4 ) + N 2 (B 4 ), appearing in (f74|), in 
terms of Ni(S 3 ) by setting 

N 2 (B 4 ) = m(S 3 ), (76) 
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and 



N (B A ) - Ni(B 4 ) + N 2 (B 4 )] = aN^S 3 ). 



(77) 



According to the above remarks, Ni(S 3 ) is asymptotically an integer multiple of 4, thus if 
we are interested to triangulations for which Ni(S 3 ) can grow arbitrarily large, it follows 
that the two parameters a and (3 necessarily are rational numbers of the form = 4 and 
a = ^ with (3 and a integers. In other words, the generic triangulations of Sf v — B A U53 
C(dB 4 ) with the joining § 3 stacked (6(3, 1) = 9/2), can be conveniently parameterized by 
setting 



and 



N 2 (B i ) ^ 
iVi(5 3 ) 4' 

N (B 4 ) - jViCgj + N 2 (B 4 ) 
iVi(S 3 ) 



a 

T 



(78) 



(79) 



where a and (3 are integers. Note that while (3 > 0, a can possibly take also negative 
values. However, if we rewrite ([T4|) in terms of such parameters 



6(4,2) = 10 



12 + 2a 



(80) 



30 + 3/3 

the kinematical bound 6(4, 2) > 4 implies 5a > 3(3, and thus a is non-negative as well. 

The parameters a and (3 so introduced are completely equivalent to the balance vari- 
ables yi and y 2 related to the cumulant action of the (fc, /) moves. Explicitly, we obtain 



%i = (^ + |/3-^Vi(S 3 )-f 
2y 2 = ( 5 -a-±P)Ni(S 3 ) + ™. 



(81) 



The Dehn-Sommerville relations for the /-vector Ni(B 4 ) allow us to express also its com- 
ponents in terms of a and (3 according to 



3N Q (B 4 ) 
3N l {B A ) 
N 2 {B A ) 
3N 3 (B 4 ) 
3N 4 {B 4 ) 



3/3 -4a 



9/3 - 10a 



N^S 3 ) + 10 
/Vi(S 3 ) + 10 



'3 + 5a 



4 

3 + 2a 



iVi(5 3 ) 



5 
2. 



(82) 
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The generic conditions 2/1 > and 2/2 > 0, (and both approaching +00), together with 
Nq(B 4 ) > 0, imply that the parameters a and 3 are related by 

\P < ol < -B, (83) 
5 4 

with (a, 3) £ N + x N + . From these remarks it follows that, as 2/1 and 2/2 go to +00, there 
are two distinct regimes for the set of triangulations considered: 

(i) If we constrain the /-vector Ni(S 3 ) of the connecting dB 4 to be O(l), then according 
to (|8l"|), a and /? go to 00 as 2/1,2/2 ~~ ► +00. From (82) we get that in this regime 



N 4 (S 4 ) ~ iV4(i? 4 ) ~ (a/6)iVi(S 3 ) 

A^(S 4 ) ^ iV 2 ( J B 4 ) = (/3/4)iV 1 (5 3 ), (84) 



where N\(S S ) is a constant. The geometrical bounds (]83| ) simply imply that as a, /3 — > 00, 
the corresponding average incidence 6(4, 2) varies between the kinematical bounds 4 < 
6(4, 2) < 5, as required. 



(ii) Conversely, if we do not constrain Ni(S 3 ) to be O(l), then according to (|8~l"D, Ni(S 3 ), 
(and hence A^3(5 3 )), is allowed to grow unboundedly large as 2/1,2/2 — ► +00. This growth, 
which corresponds to the generation of singular vertices, is possible for any finite value 
of the parameters a and 3 compatible with (jB5|). Note that if kinematically possible, 
according to (]83|), such singular triangulations entropically dominate over the regular ones 
since these latter are generated by the constrained configurations forcing Ni(S 3 ) to be 
O(l), while the former are unconstrained. More specifically, since the number of dis- 
tinct triangulations of a 3-sphere § 3 grows exponentially with A^S* 3 ), configurations with 
N^IS 3 ) as large as possible, if kinematically allowed, will dominate over configurations 
with N 3 (S 3 ) = O(l). 



The kinematical bound ( |83|) for the occurrence of singular triangulations is not trivial. 
In order to discuss its implications, let us consider the ratio between the total volume of 
the triangulated S 4 V and the volume of the ball around the singular vertex cr°, viz. 

Vol{S 4 ) Ni 12 + 2a , . 

(85) 



Vol sing (a°) N 3 (S 3 ) 9 ' 

A direct computation of the average incidence 6(4, 2), (see |80|), together with (^) immedi- 
ately shows that the smallest 6(4, 2)'s for which we may have singular triangulations occur 
for 

a = 5 + 3/i, 

3 = 8 + 5/i, (86) 

with h = 0, 1, 2, As h varies, the average incidence 6(4, 2)h and the volume ratio fl8q ) 

respectively take the values: 
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h 


P 


a 





8 


5 


1 


13 


8 


2 


18 


11 


3 


23 


14 


4 


28 


17 


5 


33 
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b(2,4) 



W ~ 4.07407 
^ ~ 4.0579 
# ~ 4.04761 
4°° ~ 4.0404 
|§ ~ 4.03508 
f§ ~ 4.03100 



Vol(S 4 ) I 



"37 

±1 

k 
1 

9 



~ 2.444 
~ 3.111 
~ 3.777 
~ 4.444 
~ 5.111 
~ 5.777 



Table 1: The smallest incidence numbers 6(4, 2)1^ and the associated singular volume 
fraction Vol(S A )/Vol s i ng (a°) as a function of the parameters a and j3. 



Vol(S 4 



Vol sing (a°) 



22 + 6/i 
9 ' 



Vol(S 4 ) 



Since the singular triangulations that entropically dominate are those for which xn~. — nrr 
is as low as possible, the smallest 6(4, 2) for which we may have generic singular triangu- 
lations with the largest Vol S i ng (o°) is 

b(2,4) sing = ^ - 4.07407..., (89) 

(corresponding to h = and y^pr^in = 22/9 ~ 2.444. . .). It is easily verified that such 
an average incidence is associated to a relative concentration of (1,5) moves versus (2,4) 
moves given by 

2/1 = 5y 2 . (90) 

It is clear from ( p7| ) that singular triangulations may appear also for smaller values 
of 6(4,2). A list of the first possible values of 6^(4,2) is provided by table |l[ These 
triangulations are less singular than the ones associated with 6(4, 2) = ^ since they 

correspond to larger values of the ratio tt-tt^i \\ , and for this reason we may be tempted 
to consider them as entropically sub dominating at least in the large- volume limit. Yet, 
this is mere apparency since their presence is particularly relevant for locating the critical 
incidence 6o and for understanding the present status of the Monte Carlo simulations. 
Moreover, as we see in the next section, these triangulations have a subtle interplay with 
the particular singular geometry dominating in the strong coupling phase of 4-D simplicial 
gravity: PL-manifolds with a single singular edge connecting two singular vertices. 

In order to get the complete geometrical picture, one has to note that for a = 2 + 8h 
and (3 = 3 + 13/i, we also get a highly degenerate configuration for which 

6 h (4,2) = ^ - 4.102564 (91) 
is a constant average incidence as h varies, whereas ) L = 16+I6fc ^ = 012,.... 
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In other words, corresponding to such value of 6(4, 2) we have distinct triangulations 
with distinct ratios t¥t~7 \\ \h but with 6(4,2) fixed. Even if this set of triangulations 

contains configurations for which y^j~^~uyU — 1-777, such a degeneration makes any 

particular configuration at fixed ) Q i entropically sub- dominating with respect to 

V Posing \& ) 

the generic configurations described by (|87|), at least as N± — > oo. 
3.4 The development of singular edges: § 4 S 

The explicit construction of the previous section may suggest that the singular trian- 
gulations we are explicitly considering are characterized by the dominance of just one 
singular vertex. Actually, as the parameters a and (3 vary, triangulations of S 4 of the 
form B 4 Ug3 C(dB 4 ) are not the only ones possible whose average incidence 6(4,2) takes 
on the value (|30|), at least as N^S 4 ) — » oo. As a matter of fact, in the infinite volume 
limit, (but not at finite volume), triangulations, with more than one singular vertex and 
with singular edges are still characterized by the average incidence (|80|). Their dominance 
in the class of triangulations considered, as N^S 4 ) — ► oo, is driven by a rather simple 
entropic mechanism which we discuss in detail in this section. 



Implicitly, the occurrence of more than one singular vertices may still be described by 
the construction B 4 U53 C(dB 4 ), since one may simply consider the new singular vertices 
and edges to be located in the ball B 4 . However, the interplay between dominance of one 
or more (edge-connected) singular vertices is most easily seen from a simple variant of the 
construction leading to (|30|). The generic singular triangulation of § 4 is still realized by 
glueing two 4- balls along an isometric § 3 boundary which is again assumed to be a stacked 
3-sphere, i.e., as B 4 S U53 B 4 . However, one of the two balls, say the one denoted by B 4 S , 
(es being an acronym for edge- singular), is no longer taken of the form of a cone C(dB 4 ) 
over the § 3 boundary, but more generally is provided by a triangulation with /-vector 





= ±J2 N z(B 3 (j)) + k 


Ni(B 4 es ) 


= ^Y, N ^B 3 (j)) + lj:N 2 (S 2 (l))+3(k 

j=l 1=1 




k k — 1 

= ^Y, N ^B 3 u)) + 2Y,ms 2 (i)) + 2(k 
6 3=1 1=1 


Ns(Bi) 


= 3j> 3 (£ 3 (j)) + |£AMS 2 (0) 
j=i z l=l 




k k—1 


N 4 (B 4 S ) 


= £iv 3 (£ 3 a)) + ]Tiv 2 (s 2 (0). 



j=l 1=1 



To grasp the geometrical origin of this /-vector imagine k distinct 3-spherical disks, B s (j), 
joined through (k — 1) § 2 -boundaries, S 2 (/); a sort of 3-dimensional peanut-shell with k- 
bulges and k — 1 necks. This gives rise to a 3-spherical peanut-shell S 3 , and we get 
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a 4-dimensional ball B^ s out of this § 3 by considering [k — l)-edges {o" 1 (0}/=i,...,fc-i n °t 
belonging to § 3 and connecting k vertices {^(j)}^!,...^ ^ §> 3 - The 4-ball B^ s is defined by 
requiring that the generic 2- spheres S 2 (l) are the links (in £> 4 S ) of the corresponding edges 
cr 1 (Z) , I = 1, . . . , k — 1. Moreover, the complex obtained by Bf s by removing the (k — 1) 
stars (in B^ s ) of the edges <r 1 (l), is assumed to be the disjoint union of k cones C(B 3 (j)) 
over the 3-spherical disks B 3 (j), with apices in the k vertices cr°(j). This construction 
can be roughly described as a 3-dimensional peanut shell containing one rather than k 
distinct 4-dimensional nuts. It is easily verified that the /- vector Ni(B^ s ) ( |92"| ) describes 
this generalized peanut triangulation and that it represents a triangulated 4-ball with 
Nq(B^ s ) vertices, k of which, {c Q (j)}j=l,...,ft are interior vertices, (i.e., o~°(j) £ dB^ s ), with 
N 3 (B 3 (j)) + N2(S 2 (j)) 4-simplices a incident on the j-th of them. The j-th of the k — 1 
interior links cr 1 (/), connects the vertex <r°(j) with <7°(j + l), and N2(S 2 (j)) 4-dimensional 
simplices cr 4 are incident on it. Thus, if some, say 1 < s < k, of the {N 3 (B 3 (j)} and the 
corresponding s — 1 of the {N2(S 2 (l))} grow with the simplicial volume of the S 4 D B^ s , 
(not necessarily with the same rate), the triangulation of B\ s just constructed contains s 
singular vertices connected by s — 1 singular edges. Note that if we take the boundary of 
this triangulated Bf s we obtain a stacked 3-sphere § 3 with /-vector 

k 

No(S 3 ) = ^N 3 (B 3 (j)) 

6 3=1 
k 

Ni(S 3 ) = ^iV 3 (S 3 (j)) 

3=1 

N 2 (S 3 ) = 2j2N 3 (B s (j)) 

3=1 
k 

N 3 (S S ) = -£N 3 (B s (j))- (93) 

3=1 

This § 3 boundary of the 4-dimensional ball B^ s may be profitably thought of as resulting 
from the connected sum, along isometric S 2 -boundaries of k distinct stacked 3-spheres § 3 , 
i = 1, ... k, to be considered as the links (in an § ) of a corresponding singular vertex. 
In this way the singular ball B^ s , (and the corresponding S 4 ), can be considered as the 
kinematical set up for discussing the interaction of k distinct singular vertices (of the 
type considered in the previous section). This picture allow also to prove an elementary 
but important result showing that, in the class of triangulations considered, the order 
of singularity of a singular edge cr 1 (/) in B^ s is subdominating with respect to the order 
of singularity of the corresponding vertices <7°(Z — 1) and cr°(Z). In other words, in the 
large volume limit, the number of 4-simplices incident on o~ l (l) grows slower than the 
number of 4-simplices incident on the vertices a (I — 1) and o~°(l). As usual, the number of 
incident 4-simplices can be considered as the (possibly) singular volume associated to the 
corresponding edge or vertex. Thus, if we denote such simplicial volumes by Fo/(o" 1 (/)) = 
#{a 4 n a l (l)}, Vol(a°(l - 1)) = #{cr 4 n <r°(Z - 1)}, and Vol(a°(l)) = #{cr 4 n a°(l)}, we 
have the following 
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Lemma 2 In the class of triangulations considered for B t 



4 

es> 



Vol((T l {l)) , , 

lim ~ ; , = 0. 94 



(Obviously the same holds with Vol(a°(l)) replaced by Vol(a°(l — 1))). In order to prove 
this result we may consider, without loss in generality, a B\ s whose S 3 -boundary consists 
of two stacked 3-spheres Si and B> 2 joined through an isometric § 2 . The more general case 
can be proved similarly without much effort. Again without loss of generality we may 
assume that the two isometric copies of S 2 along which the two S 3 and §| are glued are 
the links of a vertex in the corresponding § 3 , (as remarked in the previous paragraph, this 
can be always arranged; also it can be easily shown that there are stacked 3-spheres with 
a marked vertex whose 2- spherical link grows linearly with the simplicial volume of the 
3- sphere-see section |3.1| for an example in 4-dimensions. The above lemma states that, 
in the case of stacked 3-spheres, this large volume behavior for an § 2 cannot hold if such 
a 2-sphere is a joining neck). It immediately follows that the /-vector of dBf s = U52 S 3 , 
can be written in terms of the /-vector of §f , i = 1, 2, and of S 2 as 



iV (5 3 ) = i[JV 3 (5?) + N 3 (Sl)} - N 2 (S 2 ) + N (S 2 ) - 6 

Ni(S 3 ) = ^[N 3 (S 3 ) + N 3 (Sl)} - 4N 2 (S 2 ) + N^S 2 ) - 4 

N 2 (S 3 ) = 2[N 3 (S 3 ) + N 3 (S 3 )} - 4N 2 (S 2 ) 

N 3 (S 3 ) = N 3 (S 3 ) + N 3 (S 3 )-2N 2 (S 2 ). (95) 



Since N^S 2 ) = (3/2)N 2 (S 2 ) we immediately get 

N 3 (S 3 ) _ 3 N 3 (S 3 ) + N 3 (S 3 ) - 2N 2 (S 2 



iVi(53) 4 N^ + N.iSD-fN^)-^ 



(96) 



which implies that j^fg) > f a s long as j^gs] = 0(1) in the large volume limit (i.e., as 
N 3 (S 3 ) + N 2 (S 2 ) -> 00). Thus dB* s = U 5 2 S| can be a stacked 3-sphere if and only if 

lim = 0. (97) 

According to the above remarks N 2 (S 2 ) = Vol^ 1 ^)) and N 3 (S 3 ) + N 2 (S 2 ) = Vol(a°(i)), 
and we can write (|97|) as 

, , V t iaHt) \ f7T\ = 0, (98) 

from which the lemma follows. 



This latter result only implies that the singular volume of the edge cannot grow linearly 
with the total volume of the ball B^ s (and of the resulting § 4 -see below). As we have seen 
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in the previous section, a linear growth is instead typical for the singular volume associated 
to the vertices. It should be stressed that a subdominant rate of growth, (say with some 
fractional power of the total volume of B^ s ), is well in agreement with (p7| ). As a matter 
of fact subdominant powers for the volume growth associated with a singular edge are the 
ones typically experienced in numerical simulations ||. 

Note that triangulations with Vol(a l (i)) as large as kinematically possible, thus grow- 
ing with [N A {B^ S )\ 5 for some < 5 < 1, entropically dominate over triangulations of 
Bg S with Volfa 1 ^)) = 0(1). This remark follows as a direct consequence of the fact 
that triangulating B^ s under the hypothesis N<2(S 2 ) = 0(1), while sufficient to assure the 
validity of (|97|), it is not a necessary condition. It generates a subclass of constrained 
configurations in the class of triangulations of B^ s considered. Conversely, triangulations 
with N2(S 2 ) oc [A?4 (£?£,,)] 5 are, according to (|97|), unconstrained, and as such much more 
numerous at least in the large volume limit. 

As in the previous section, we obtain a 4-sphere §g S , (es again for edge- singular), by 
glueing a generic triangulated ball B 4 with stacked S 3 boundary to the singular B 4 S defined 
by (^), viz., §f s ~ B 4 U S 3 B 4 S . It is easily checked that the /-vector of such an S 4 S is 
given by 

N (S 4 S ) = N (S 3 ) + N (B 4 ) + k 

miSi) = N 1 (S 3 ) + N (S 3 )+N 1 (B 4 ) + W N ^ s2 ( l ))+ 3 ( k - 1 ) 

1 i=i 

k-l 

N 2 (S A es ) = N 2 (S 3 ) + N 1 (S 3 )+N 2 (B 4 ) + 2j2N 2 (S 2 (l)) + 2(k-l) 

i=i 

mst) = N 3 (S 3 ) + N 2 (S 3 ) + N 3 (B 4 ) + IJ2n 2 (S 2 (1)) 

/ 1=1 

k—1 

mS 4 es ) = N 3 (S 3 ) + N^B 4 ) + J2N2(S 2 (l)), (99) 

i=i 

where Ni(S 3 ) denotes the /-vector of the joining stacked 3- sphere d(B 4 s ) ~ S 3 , and Ni(B A ) 
is the /-vector of the interior of B A . According to (|97|), N 2 (S 2 ) / N 3 (S 3 ) is asymptotically 
o(l), thus, in the large volume limit, the average incidence 6(4,2) of such a triangulated 
§g\, is still provided by the expression ( |80[ ) introduced in the previous section, viz., 

^U 6 < 4 ' 2 >'* = 10 'iw (100) 

where the two parameters (3 and a are again defined by (78) and (|79|), respectively. Before 



we proceed any further, we should emphasize that ( |100| ) strictly speaking only holds in the 
limit iV^S" 4 ) — > oo, and that at finite (but large) volume N^^S^), we have 6(4, 2)| S 4 > 
6(4, 2) with 

,/ i 12 + 2a ZZ^'I N 2 (S 2 (l)) , . 

t ' 4 - 2 >^= 10 -30T3g + " N 3 (S1) • < 101 » 
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for a suitable a- and /3-dependent constant n > which can be easily worked out. For 
instance, for the relevant case k = 2, (i.e., two singular vertex connected by a sub-singular 
edge), we get to leading order in N2(S 2 )/N 5 (S 3 ), 



6(4,2)| 



10 



12 + 2a 
30 + 3/3 



+ 10 



6 + 3/3 - 4a 
100 + /3 2 + 20/3 



jV 2 (S 2 ) 
iV 3 (5 3 ) 



(102) 



Since according to lemma |2] the ratio ^=j^rg£j y °" can g° to zero, in the large volume 
limit, as slowly as A^S 3 )* 5-1 for some < 5 < 1, we get the 



the singular-vertex triangulations 



Lemma 3 At finite volume N4(S 
C(dB 4 ), considered in section 3.S, are closer to the kinematical boundary 6(4,2 



B s 



U 5 3 

= 4 



than the edge-singular triangulations § 4 S ~ B U53 B\ 



We stress that this result does not imply that the singular-vertex triangulations ~ 
B A U53 C(dB A ) entropically dominate in the large volume limit. For, according to ( |100| ), 
the edge-singular triangulations become more and more important as the volume increases, 
and eventually in the infinite volume limit the triangulated spheres S 4 S enter in full entropic 
competition with the triangulated Sf„ considered in the previous section. Actually this 
entropic competition comes into play quite rapidly as the volume increases. For instance, 
from ( 102 ), one gets that, for the dominating configurations at h = 0, 



... 110 100 

6(4,2) ,4 ~ 

v i\b M 27 324 



N 2 (S 2 ) 



(103) 



Numerical simulations at A^(5 4 ) = 32000, (see e.g., || and show evidence that 
N 2 (S 2 )/N 3 (S 3 ) < 1/10, thus the average incidence 6(4, 2)1^ of Sj s differs (at h = 0) from 
the average incidence 6(4, 2)| 5 4 s of Sf, s by less than 3/100. Therefore it is important to 
understand how, as -/V^S 4 ) increases, the k distinct singular vertices (and the correspond- 
ing k — 1 subsingular connecting edges) interact among them, and which configuration 
actually dominates in the large volume limit. 



As we have seen in section |3.3| , the various singular triangulations of the 4-sphere 
considered there are parameterized by the ratio between the total simplicial volume of the 
given Sg V and the simplicial volume of its singular part, (see fl85|) and (|8q) ). If we consider 
a similar ratio also for S 4 S , i.e., 

Vol{Sj s ) = N 4 (Sj) 

Vol(sing) iV 4 (B 4 s ) ! 1 ' 

then, as is easily verified, this ratio is still provided, in the large volume limit, by (|85|). It 
follows that the entropic comparison between the single singular vertex triangulations § 4 „ 
and the multiple singular vertices triangulations § 4 S should be carried out at a fixed value 
of the ratio Vol(S^ s ) /Vol(sing) = const. = Vol(Sg V )/Vol S i ng (a°). 
In our case 

k k—1 

N,{Bi) = $>3(B 3 (j)) + £^(S 2 (Z)) 
i=i (=i 
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fc-1 



iV 3 (5 3 ) + ^iV 2 (5 2 (0). (105) 



1=1 

According to the remarks following lemma ||, unconstrained triangulations of Sg S generally 
have J2ti N 2 (S 2 (l)) = 0{N 3 (S 3 ) 5 ) for some < 6 < 1. Thus 

lim N,(Bi)/N 3 (S 3 ) = 1, (106) 

AT 4 (S 4 )^oo 

and working at constant ratio ( |104| ), (in the infinite volume limit Na(S^. s ) — > oo), implies 
that we have to consider triangulations of -B 4 ,, with 

AT 3 (S 3 ) = Aj • A^(S 4 ) (107) 

and 

fc-i 

A 2 < ]T N 2( s2 ( 1 )) < A 3 ■ N 3 (S 3 ) 6 , (108) 
i=i 



for some positive constants A\, A%, and A 3 . 

Guided by these considerations we can easily get a set of entropic rules for deter- 
mining which configuration dominates in the set of singular triangulations of S 4 . We 
start by an obvious adaptation of an argument in 0, according to which the number 
of distinct triangulations associated with a singular vertex, (the local entropy of the ver- 
tex), is provided by the number of distinct triangulations of the link of the given vertex. 
The link, link(o~°(j)), around the j-th singular vertex o~°(J) £ B^ s , is a 3-sphere § 3 (j), 
and any two such links, S 3 (j) and S 3 (j + 1), associated with two singular vertex con- 
nected by a singular edge o~ l {j), have a non-empty intersection § 2 (j), (the link of the 
connecting edge Thus, the inclusion-exclusion principle implies that the number, 

Card[B* s (S 3 (l), S 3 (k); S 2 (l), ...,S 2 (k- 1))], of distinct triangulations of B^ s with 
given singular vertices {5 3 C7)}i=i,...,fc and given singular edges {<S' 2 (Z)}z=l,...,fc-l is pro- 
vided by 

a , <. 9 nf-i Card[S 3 (j)} 

Card[Bi(S 3 (l), S 3 (k); S 2 (l), S 2 (k - 1)))] = ^J^j , (109) 

where Card[S 3 (j)\ and Card[S 2 {l)] respectively denote the number of distinct triangula- 
tions of the 3-spherical links of the j-th. singular vertex and of the 2-spherical singular link 
of the Z-th singular edge. Since each S 3 (j) is a stacked 3-sphere, (hence with an average 
incidence 6(3, 1) = |), the microcanonical partition function (||) immediately provides the 
leading order asymptotics both for Card[S 3 (j)] and Card[S 2 (l)], viz., 



Card[S 3 {j)] N3{sS{j)) 



»i — 



(i( 3,i)_g + i)K3,iH + ii w '( s3 ^ 



{b{3, 1) - g) fe (3,lM 

where b(3, 1) = §, q = 3. Since A^i(5 3 (j)) = ^N 3 (S 3 (j)), we get 

-(5)5/21 l^3(S 3 (i)) 



(110) 



( 3 )3 /2 



(111) 
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CardlS 2 



\N 2 {S 2 (l))»l 



(6(2,1) -g+ 1)6(2,1)-?+! 



= 2 + lN 2 (S 2 (l)) 


(H|) provides 




" 4 4- 


i7V 2 (S 2 (/)) 




3 3 


• (112) 



(6(2,1) -g) fe (2,l)-g 

Thus, by setting (7(2) = [4V3 3 ] 1 / 2 and C(3) = [(5/2) 5 / 2 /(3/2) 3 / 2 ] 4 / 3 , we eventually get 



Card[B* a (S 6 (l), S 3 (k); S z (l), ...,S 2 (k- 1))] ~ 

k \ /k—1 \ 

N 3 (S 3 (j)) In C(3) - f £ iV 2 (5 2 (0) ) InC(2) 



cxp 



(113) 



Since 



£ iV 3 (5 3 (i)) = iV 3 (S 3 ) + 2J2 N 2 (S 2 (l)), 
j=i 1=1 



(114) 



where S 3 = <9-B 4 s ^ s l ne stacked boundary of U 4 S , we can rewrite (113) as 

C(3) 2 



Card[Bi(S 3 (l), . . . ; S 2 (l), ...,)]- C(3)^ 3 ) • 



C(2) 



E*ri ljv 2( s2 (0) 



(115) 



(by exploiting ( |105| ) this expression can be also rewritten in terms of N^(B^ S )). Since 
C(3)/C(2) > 1, we have that triangulations of B\ s with large Y$=i N 2(S 2 (l)) are dominant 
in the infinite volume limit. This implies that the simplicial volume of the k — 1 edges 
connecting the k vertices is as large as possible. Note that ( |115| ) does not depend on 
the particular S 3 (j) or S 2 {1) but only on the fixed quantities N 3 (S 3 ) and ^=1 N 2(S 2 (l)) 
determining the ratio between N^S 4 ) and the volume of the singular part B^ s of § 4 , (see 
( jLQjfl and (|105D). 

Thus, among all possible triangulations with k distinct singular vertices connected by 
k — 1 distinct edges, those entropically favored, as k varies, are the less constrained ones, 
namely triangulations with just one singular edge connecting two singular vertices: the 
triangulations of B^ s with k = 2. For such triangulations the S 3 links of the singular 
vertices and the S* 2 link of the connecting edge are as large as kinematically possible. 



Note that for the triangulated B 4 = C(dB 4 ) considered in section 3^ we have 

Card[C(dB 4 )] ~ C(3) N ^ s3 \ 
and in the large volume limit 

Card[Bi(S 3 (l),...;S 2 (l),...,)]^ 



(116) 



C(3) 



N 3 (S 3 ) . 



g(3) 2 
C(2) 



Efj> 2 (s 2 (o) 



> 



Card[C(dB 4 )} ~ C^)^ 53 ). (117) 



Since, as A^(5 4 ) increases, the triangulations S 4 S enter more and more in entropic compe- 
tition with the single singular vertex triangulations § 4 „, ( |117| ) directly implies the following 
basic result 
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Lemma 4 For a given ratio 

Vol(St) _ Ni(Si) _ 22 + 6/j 



(118) 



Vol(sing) N 4 (Bi s ) 9 ' 

with h = 0, 1, 2, . . . ,, the singular triangulations of § 4 which are closer to the kinematical 
boundary 6(4, 2) = 4, and which entropically dominate in the large volume limit N^S 4 ) — > 
oo ; are realized by triangulations with one sub-singular edge connecting two singular 
vertices, and are characterized by the average incidence 



The last part of this lemma, concerning the /i-parameterization of the singular trian- 



gulations, is an immediate consequence of the expressions ( |100| ) and (101) for the average 
incidence of §f s and of the results of section [T^. Results which characterize the sets of 
value of a and (3 giving the closest approach of 6(4, 2) = 10 • 39+3^ to the kinematical 

boundary 6(4, 2) = 4 as the ratio y^^g) varies. 

The geometrical analysis just discussed and the Lemma ||| appear in good qualitative 
agreement with the picture which emerges from recent Monte Carlo simulations |7|] con- 
cerning the study of singular structures in 4D simplicial gravity. According to such a 
numerical analysis there are, at finite volume, two pseudo-critical couplings (and hence 
corresponding pseudo-critical incidences 6(4, 2)) separately associated with the creation of 
singular edges and singular vertices. This behavior seem to correspond to the different 
entropic relevance of the single singular vertex triangulations Sj v and of the singular edge 
triangulations Sf s discussed above. In the simulations the two pseudo-critical couplings 
lock into a single critical point in the large volume limit. This merging appears to be re- 
lated to the full entropic competition between S^ v and §g S which dominates our geometrical 



picture in the infinite volume limit. Explicitly, the average incidence 6(4, 2) 1 54 , (see 101), 
is slightly larger (at finite volume) than 6(4, 2)|g4 . Thus, if we apply formula ( jl4|) relating 
the average incidence 6(4, 2) to a value of the coupling &2, we find that the set of ^(S^J's 
corresponding to 6(4, 2)|g4 , (as h varies), is slightly smaller than the corresponding set 
of /^(S^J's associated with 6(4, 2)\ S 4 . Anticipating the analysis of section |], this remark 
implies that there are indeed two pseudo-critical points respectively associated with edge- 
singular §g S and vertex-singular S^ s triangulations, say kSf 1 * (S%g] N&) and kg" 1 * (Sy S ; N4) , 
with 

kT\Si s ;N A )<kT\St s ;N A ), (120) 

and coalescing in just one critical point as N4 gets larger and larger. Obviously, what 
one actually sees at a given finite volume mostly depends on the rate A*2 (S^/iV^S 13 ), (see 



101), which controls how fast the two average incidences 6(4, 2)|g4 and 6(4, 2)\ S 4 approach 
each other. On this rate we are not yet able to say anything substantial. As recalled, 
(see (|103D ), computer simulations indicates that at relatively large volumes, (tipycally 
N 4 = 32000), the term N 2 (S 2 )/N 3 (S 3 ) is already so small that 6(4,2)1^ ~ 6(4, 2)1^ 
up to a few percent, and edge-singular triangulations are to all effects as close to the 
kinematical boundary 6(4, 2) = 4 as the are. Thus they do entropically dominate. 
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3.5 The characterization of the critical incidence 

Since in the infinite volume limit both singular configurations S 4 V and § 4 S are characterized 
by the same average incidence (]80|), we can use indifferently both for characterizing the 
critical incidence 60 (4) signaling the closest approach of generic singular triangulations to 
the kinematical boundary 6(4, 2) = 4. The single singular vertex configurations Sf v are 
somehow easier to handle than S 4 ,,, thus for definiteness we describe the characterization 
of the critical incidence (and the corresponding critical gravitational coupling) by referring 
explicitly to S 4 „ ~ B U53 C(dB 4 ). In any case, one should keep in mind that the extension 
of the analysis to § 4 S can be carried out without difficulty along the same lines. 

How can we characterize the critical incidence 60 (4)? A glance at table [I] clearly shows 
that, as v V i° l ^ S ( \\ increases, the values of 6(4, 2)L are very close to each other. This remark 
implies that triangulations with 6(4,2)1^=0 = 110/27, even if entropically dominating in 
Sg V ~ B 4 U53 C(dB 4 ), cannot be taken as the mark of the real critical incidence. As a 
matter of fact, for values of h close to the leading configuration at h = 0, there can be 
statistical competition between such singular triangulations, at least as N4 — > 00. The 
critical incidence 60 is actually obtained by averaging the distinct 6(4, 2) | ^ 's over the set 
of corresponding singular triangulations. 

To characterize such average we exploit the fact that the singular triangulations we are 
considering have their singular part constructed as a cone over a stacked 3-sphere § 3 . If we 
join, through the identification of a marked <r 3 £ S 3 , two stacked 3-spheres, § 3 (1) and § 3 (2) 
we get another stacked 3-sphere S 3 (3) = § 3 (1)#§ 3 (2), and all (voluminous) stacked spheres 
can be obtained in this way. Thus, if we construct the cone over this connected sum of 
stacked 3-spheres we can sweep all possible voluminous (i.e., large N4) singular triangula- 
tions of the type we are considering. Explicitly, let us denote the singular triangulations of 
S 4 , obtained from the stacked 3-spheres S 3 (l) and S 3 (2), by S 4 (l) = B 4 (l) U s z {1) C(§ 3 (1)) 
and S 4 (2) = B 4 {2) U 5 3 (2) C(§ 3 (2)), respectively. If S 3 (3) = § 3 (l)#/§ 3 (2), where / is an 
homeomorphism between two marked cr 3 (l) £ § 3 (1) and <r 3 (2) £ § 3 (2), then 

§ 4 (3) = S 4 (l)# / ,§ 4 (2) = (B\l) U f B\2)) U 5 3 (3) C(S 3 (3)), (121) 

where /* is the extension of / to the cone over the marked a 3 , and every singular trian- 
gulations of § 4 over a stacked 3-sphere can be obtained in this way. 



The analytical counterpart of ( |l2l| ) follows directly from the last of relations (82) 
characterizing the /-vector of the ball B 4 as the parameters a and /3, (thus h), vary. From 
it we get 



N,[B 4 (l)U f B 4 {2)) 



13 + 6/i 



9 



A^i(5 3 (l)) + 



13 + 6/i 



A^i(5 3 (2)) 



N 4 [B 4 (1)] + N A [B%2)], (122) 



where we have discarded costant terms which are o(l) in the large N4 limit. To exploit 
this information let 



T h [Vol(B 4 ) =N} = Card { S 4 : V f^ S ^ = Q l±^ . V ol{B 4 ) =N\, (123) 

I Volsing^) 9 I 
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be the cardinality of the set of distinct singular triangulations of the ball B 4 , constructed 



over a stacked § 3 , with given ratio 



Vol(S 4 ) 
Vol s i n g{a u ) 



and N 4 (B 4 ) = N. According to the behavior 



of this set of triangulations under the connected sum we have 



T h 



Vol{B 4 ) = JV(1) + JV(2) 



Vol(B 4 ) = N(l] 



T h 



Vol(B 4 ) = JV(2)1 . (124) 



It is easily verified that this relation implies that the leading asymptotics of Th(Vol(B 4 )) 
is provided by 



T h (Vol(B 4 )) = c(B 4 ;h) 



Af 4 (B 4 ) 



(125) 



where In c(B ; h) is the specific entropy for the generic a € S 4 S . 

Since there is a unique cone C(§ 3 ) over the stacked sphere boundary dB 4 ~ S 3 , (|l2 
provides also the leading exponential asymptotics to the number of distinct triangulations 
of § 4 S with given N 4 and given h, viz., 



Card{si s } oc c(B 4 ; h) 



(126) 



Actually, when h » 1 and N^^S 3 ) = 0(1), for each triangulation of S 3 , there can be a 
worth of Aut(S 3 ) inequivalents cones, Aut(S 3 ) denoting the automorphisms group of the 
given triangulation, for simplicity we disregard here these correction factors. Note also that 
the above construction applies to the edge-singular spheres S 4 S with minor modifications. 
According to ( 34 ) , 

N 4 (B 4 ) = N 4 (S 4 ) - N 3 (S 3 ) = N 4 ( 

thus we get that to leading order 

Card{S 4 vs }=c(B 4 ;h) N ^- N ^ 



13 + 6/i 
22 + 6/i : 



(127) 



s(h) N ^ s4 \ 

where we have introduced the specific entropy, In s(h), of a a 4 £ §' 

In Car d{S 4 s } 



(128) 



according to 



In s{h) 



lim 

iV 4 (S 4 )^oo 

13 + 6/i 
22 + 6/i 



N 4 (S 4 ) 
\nc{B 4 -h). 



(129) 



In order to characterize lns(/i), note that triangulations of the form S 4 S describe, 
for h = 0, the generic singular triangulations of § 4 realizing the closest approach to the 
kinematical boundary 6(4, 2) = 4. Conversely, and as already stressed, the triangulations 
S 4 S reduce, as h — > oo, to the generic (branched polymer) triangulations of S 4 , (with a 
rooted a 4 ). These remarks imply that corresponding to h = and h = h max we must have 



In s(h = 0) 
In s(h = h max ) 



lnc(S 4 ;h = 0) 
In c(S 4 ; h = h m 



(130) 
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where h max is characterized by the value of the ratio (|8^) evaluated for the smallest 
possible iV^S" 3 ) = 5, i.e., h max = ^j-/V 4 — and where lnc(5 4 ;/i) is the specific entropy 
associated with the microcanonical partition function (|lO|) , i.e., 



c(S 4 ;h) 



(K 4,2)-2)W- 2 l 10/6(4 ' 2) 



(131) 



(6(4,2)-3) 6 (4,2)-3 

with 6(4,2) = 10 • (22 + 6/i)/(54 + 15h), (the actual specific entropy contains a constant 
factor which is of no relevance for the present considerations-see (10)). 

Since c(S 4 ;h) is a slowly varying function of h, the specific entropy lns(/i) can be 
characterized as the convex combination of In s(h = 0) and In s(h = h max ) over the interval 
< h < h max , viz., 



In s{h) 



h 



hr. 



■ In s(h = h Ti 



+ 1 



h 



h, 



-) \ns(h = 0). 
/ 



(132) 



In other words, we are considering lns(h) as the convex combination of the extreme 
pure phases (h = 0: crumpling, and h — > oo: branched polymer). A straightforward 
computation provides 



s{h) = c(S 4 ;h = 0) 



c{S i ;h = 0) 
c(S ; h = h max ) 



3AT 4 



Since in the large N 4 {S A ) limit, ln[c(5 4 ; /i = 0)/c(5 4 ; /i = h T 
for the leading asymptotics 



Card{3t s } = c(S A ; h = 0) Ni 



e s . 



(133) 

0.06 we eventually get 
(134) 



It is worth stressing that a completely analogous result holds for Card{§f s }, since, as 
N4 — ► 00, the set of edge-singular triangulations, (with one edge connecting two singular 
vertices), S 4 s |fc=2, is as close to the kinematical boundary 6(4, 2) = 4 as the triangulations 
S vs . The two class § 4 S and § 4 S only differ in the subleading asymptotics. 



According to (p.34|) , the average value of 6(4, 2)\h over the set of singular triangulations 
considered is given, in the large N4 limit, by 



(bU2) )\ n=oKi,2)\ h exp[ 

2^h=o eyL Pi 5 J 
By approximating the numerator with an integral, we get 



(135) 



{b(4,2) sing )\ hn 



4 + — 
15 



4 e-^E^-E^ 



+ 



18 - 
25' 



5(1 



(136) 



where E\{x) is the exponential integral function. In the large volume limit h r , 
and the above expression reduces to 



00, 



(6(4,2; 



sing I 



4 « ,18, 
75 U 25' 



4.0394361235. 



(137) 
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As stressed, a similar analysis carried out for the class of singular triangulations S 4 ,, would 
provide the same (6(4, 2) s i ng ). It follows that, as A^S 4 ) — > oo, (137|) is the value of 



the incidence 6(4,2) statistically dominating in both sets Sj. v and S 4 S . As argued in the 
previous sections, these triangulations are the ones characterizing the smallest possible 
6(4, 2) marking the onset of the dominance of singular geometries. Thus, we can identify 
(6(4, 2) s i ng ) with the critical incidence 6o (see section |2.1| ) characterizing the transition 
between the weak and the strong coupling phase of the theory, i.e., 

6 (4) = (6(4, 2) smg ) ~ 4.0394361235. (138) 



Together with the critical incidence (6(4,2)) it is worthwhile to compute the infinite 
volume average, over the set of singular triangulations Sf v or S 4 ,,, of the local volume of 
the singular part of the triangulation, Vol(sing). Note that for the class of triangulations 
S 4 „, Vol(sing) = Vol (a ), whereas for the triangulations of § 4 S dominating in the infinite 
volume limit, we have 

Vol(sing) ~ 2Vol(a°), (139) 



since according to the remarks of section |3.4| and lemma in such a limit, triangulations 
with just two singular vertices (connected by a sub- singular edge) dominate. 

For both class of triangulations Vol(S^) /Vol(sing) = (22 + 6/i)/9, and the required 
average is provided by 

Vol(sing) E/gf ex P[-f]6feT22 mm 

\ vdm Jlfc- £fc-exp[-J] ' 1 ' 

(Strictly speaking, this ensemble average explicitly refers to the single singular vertex 
triangulations S 4 ^, however, as stressed before, this ensemble average differs from the § 4 S 
ensemble average by corrections which vanish as N^S 4 ) — ► oo). 

By approximating as usual the summatories with an integral we get 

Voljsing) 3e n / 15 f ,11 „ ( h max 11 

( voi(s^) )lh — " io(i - e -^/5) • ^TE' ~ E ^— + TE )] - (141) 

According to ( |139| ), we get for the average local volume of the (two) most singular 
vertices, the explicit expression 

which, in the infinite volume limit, reduces to 

Q p H/15 11 

(Vol(a )) = ^-Eii-) ■ N 4 . (143) 

Note that the value of the critical average incidence (6(4,2)) ~ 4.03943... shows 
that the leading configurations contributing to the singular geometry of S 4 S are, loosely 
speaking, those for which h < 6, (see table |2|). Thus, a rough indicator of what is the 
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average singular volume for 5(4,2) sufficiently smaller than (6(4,2)) ~ 4.03943..., (viz., 
when in the polymeric phase), can be obtained by considering the average 

j Voljsing) ^ _ Efe|tT ex P[~f ] 6FF22 nAA \ 



Explicitly we get 



o 29/15 OQ 

(Vol(a°)) poly = • JV4, (145) 



which can be interpreted as the contribution to (Vol(a )) coming from the non-singular 
geometries in §g S . 



4 The critical coupling k^ 11 

The kinematical picture which emerges from the above analysis is immediately connected 
to the thermodynamical behavior of 4D- dynamical triangulations by recalling the results 
of section according to which, as k 2 varies the distribution of triangulated manifolds is 
strongly peaked around triangulations with an average incidence given by 3[A(k2)/ (A(k2) — 



1)], (see (14)). Thus by solving for k 2 the equation 

(b(4,2U g )=3( (146) 

we get an estimate of the value of k<i corresponding to which singular triangulations start 
dominating the canonical partition function (Q) in the infinite volume limit. Recall that 
singular triangulations are those characterizing the sub-exponential sub-leading asymp- 
totics (see Th.5.2.1, pp.106-118 of |§) 

W[N 2 ,b(4,2)] - (147) 



>4&(4,2))AT : 



{b-q + l) b ^ +1 



ib-q) 



b-q 



N 2 



with m(6(4, 2) > 0, (see @ for the general expression; the above expression can be 
obtained from (]9|) by setting n = 4, 02 = 0, r(6) = 0, and D = since we are considering 
§ 4 , we have also dropped a few inessential constant terms). Thus we can identify the k 2 
solution of equation ( |146j ) with the critical value, k^, of the inverse gravitational coupling 
marking the transition between the strong and weak coupling in 4D-simplicial quantum 
gravity. 



Introducing in ( |146[ ) the values {b(4,2) s i ng ) ~ 4.0394361235 obtained above for the 
kinematical bound controlling the occurrence of generic singular triangulations, we get for 
the critical coupling the explicit value 

fcf u ~ 1.3093. (148) 
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h 


b(2,4) 


Vol(S*) | 
Vo/ s i„<,(o- u ) 1" 


ju/i 





W ~ 4.07407 


2.444 


~ 1.24465 


1 


W ~ 4.0579 


3.111 


~ 1.2744 


2 


# ~ 4.04761 


3.777 


~ 1.2938 


3 


1? ~ 4.0404 


4.444 


~ 1.30746 


4 


fff ~ 4.03508 


5.111 


~ 1.31762 


5 


§§ ~ 4.03100 


5.777 


~ 1.32545 



Table 2: Some of the values of k\ obtained by solving equation (|146|) for 6^(4,2) as h 
varies. Such values appear strikingly near to the values of the pseudo-critical points found 
in Monte Carlo simulations as the size of the triangulations considered is increased. 

4.1 A model for pseudo-criticality at finite N^S 4 ) 

It is very interesting to compare the value for k^, already in very good agreement with 
what is found by means of Monte Carlo simulations, with the other fc^'s obtained by 
solving equation (|146 ) with the left member (6 s j ng (4, 2)) replaced by the values 6/j(4, 2) 



provided by (^7j). In this way we get table || 

According to the remarks in the previous paragraph, k%, h = 1,2,..., can be interpreted 
as the values of the inverse gravitational coupling corresponding to which the sub-leading 
singular configurations comes into play. In other words, corresponding to such values of 
hi there are distinct peaks in the distribution of singular triangulations of § 4 S . The leading 
peakis at &2 = &2™* — 1-24465, this corresponds to the dominance of singular triangulations 
for which y~f~f~n^ \h = 22/9; the first sub- leading peak occurs at &2 = k\ =1 ~ 1.2744, 

corresponding to the sub-dominance of singular triangulations for which v V , ol ^ S , \ \ \h = 
28/9; the second sub- leading peak occurs at = &2 =2 — 1-2938 and is associated with the 
sub-dominance of singular triangulations for which y^°, \h = ^^/^i an d so on. In the 
large ./V4 limit there is enough phase space for having all such peaks contributing to the 
partition function of the theory, and the presence of the sub-dominating peaks lowers the 
critical incidence from its bare value 6(4, 2)|^ = o to (6(4, 2) s j ns ), and shifts the critical k^ 11 
from the bare value 1.24465 to its effective value k^ 11 — 1.3093. Using a field-theoretic 
image, one may say that in the large volume limit the fluctuations associated with the 
various sub-dominating peaks in the distribution of singular triangulations dress the bare 
critical incidence. 

Conversely, at a finite value of N4 one expects that the resulting average (6(4, 2) S i ng ) , 
computed from ( |135| ) with h < h(N^) < h max , for some h(N^), is larger than the limit- 
ing value (6(4, 2) S i ng ). Corresponding to this (6(4, 2) S i ng )(N±) one gets an JV4- dependent 
pseudo-critical point k^iN^) smaller than the actual k^ 11 ■ Roughly speaking, at finite 
volume, there is no phase space available for having all sub-dominating peaks competing 
with each other according to their relative entropic relevance. Moreover, at finite volume 
we should distinguish which kind of singular geometry we are dealing with. According 



to lemma || and (101), the average incidence is larger for the edge-singular triangulations 



5g S than for the single singular vertex triangulations § 4 . Thus, corresponding to S 4 or S 4 S 
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we should get a slightly different sequence of pseudo-critical points, (according to ( |146| ), 
k^' 11 {N ^)\{S 4 S ) < kf^iN^KS 4 )). A difference which fades away as the volume increases. 



In order to make contact with numerical simulations is worthwhile to develop an ana- 
lytical model taking care of these finite size effects. Again for simplicity, let us limit our 
analysis to the vertex singular triangulations S 4 S , with the understanding that what we say 
can be easily extended to the edge-singular triangulations S 4 S with minor modifications. 



The starting point of our analysis is the entropic formula (134) expressing, as h varies, the 
entropy of the triangulations S 4 S as convex combinations of its extreme two pure phases 
associated with crumpling {h = 0) and polymerization (h = h max — > oo). Rather than use 
directly (|134[) we should refer to the conditional entropy 



CardS 4 s 
CardS 4 



(149) 



which provides the contribution of the triangulations S 4 S to the set of all possible trian- 
gulations of § 4 , at fixed volume. 

From (134) and (S) we get, to leading order in the large N^(S 4 ) limit, 



CardSt 



CardS 4 

where f2 is the /i-dependent constant 
. c{S 4 ;h = 0) 



n 



c{S 4 ;h) 



33.97082 



(6(4,2) -2)^M ] - 10/h(4 ' 2) 
(6(4,2) -3) & ( 4 . 2 M 



(150) 



(151) 



with 6(4, 2) = 10 • (22 + 6/i)/(54 + 15/i). 

The expression ( |150| ) for the conditional entropy, holds at finite, sufficiently large 



CardW 



< 1, it implies that, at finite volume, triangulations § with 



N±(S ), and, since 

h » N^(S 4 ) In Q are entropically suppressed. This remark implies that the configurations 
S 4 S , which actually contribute in characterizing the critical incidence, have an entropic cut 
at some value of h, say h(Ni) = 0(Ni(S 4 ) In Q). The specific entropy In c(S 4 ;h) of {§ 4 } 
changes very slowly with h, thus at finite N^S 4 ) = N, we may tentatively write 



CardS 4 s 
CardS 4 



"o 

N 



N ^ (152) 



for < h < h{N) whereas 



CardS 4 s \ _ Ni 



CardS 4 



^h-h > ( 153 ) 

fl — Umax K ' 



Lmax 

N 



for h(N) < h < h max , and where f2o = 1 + e > e > 0, is a suitable constant not differing 
much from 1, (according to (pip, Sl\(h = 1) ~ 1.01234, and Q\(h = 10 5 ) ~ 1.0615). In 
other words, we are assuming that for < h < h(N) the system may exist as a mixture of 
its two extreme pure phases, whereas for h > h(N) it collapses into its branched polymer 
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phase. It is worthwhile stressing that more realistically one may consider, in place of 
(152), a convex combination of the extreme phase h = and the (non extreme) phase 



corresponding to h = h(N). By exploitying ( [L33D , this prescription can be worked out 
without difficulty, however it gives rise to a rather complex scaling behavior of the resulting 
entropy. Moreover, the fact that c(S 4 ; h) is a slow varying function of h, makes, as we 
shall see, the simpler ( |152| ) quite accurate and much easier to handle. 

A qualitative characterization of h(N) as TV varies can be easily obtained by the obvious 
scaling properties of ( |152j ). If we consider triangulations § 4 S with two distinct volumes, 
say N 4 (S 4 ) = N(l) and N 4 (S 4 ) = N(2), then 



' CardSt* \ f CardSt 



(154) 



\ CardS 4 I r „,„, I CardE, 4 I T „,„, 

\ / N 4 =N(1) \ / N4,=N(2) 

provided that h(N) scales with N according to 

h(N(2)) = h(N(l)) + 5[iV(2) - JV(1)] In fi„- (155) 

This scaling relation implies that h(N) has a linear dependence on N^S 4 ) according to 

h{N A ) = 5N 4 lnn (156) 

where £ is a suitable constant. This rather simple argument does not yet provide the 
actual value of the constants f^o and £, however confrontation with numerical data at 
A^4(5* 4 ) = 32000 indicates as reliable candidates the values 

5 In r?o : 



16000 

£ = -1. (157) 

Note that the above condition for Qq implies Qq ~ 1.0000125, a value which is perfectly 
consistent with the above characterization of f2, (see ( |151|) . It also indicates that the 
triangulations § 4 S , (actually the entropically dominating § 4 S ), do saturate the possible set 
of triangulations of § 4 in the strong coupling phase. 

The iVVdependent pseudo-critical incidence (6(4, 2) s i ng ){N 4 ) and the associated pseudo- 
critical point kf^iN^) Ccin be Gcisily obtained from ( ]13€| ) by replacing h max with h(N), 
viz., 

15 5(1 -e—T 1 ) 

\/AT.\ !v, „f /UtA OA . \ 

sing I • 



and by solving for k,2 the equation ( |146|) with (6(4, 2) s i ng ){N 4 ) in place of (6(4,2] 



By exploiting these results we get an overall analytic picture of the large volume 
behavior of 4-dimensional simplicial quantum gravity which is in a surprising agreement 
with the Monte Carlo simulations of the real system 0] . 
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5 Comparison with Numerical Work 



At this stage it is indeed useful to discuss the status of our geometrical results in the light 
of the most recent numerical work. This comparison is particularly important since , as 
recalled in the introductory remarks, the current perspective on 4-dimensional simplicial 
quantum gravity has undergone a rather drastic change. As a matter of fact, recent Monte 
Carlo simulations seem to accumulate more and more evidence for a first order nature of 
the transition separating the strong and the weak coupling regime of the theory. Taken 
at face value this result suggests that dynamical triangulations is not likely to be a viable 
model of quantum gravity unless one adds additional terms to the action. It is perhaps fair 
to say that the geometrical analysis of the previous paragraphs bears relevance to such an 
issue. The characterization of the critical coupling k^ 11 and the existence of entropically 
sub- dominating peaks in the distribution of singular triangulations strongly indicates that 
this geometrical picture may be responsible for the phenomenology we see in numerical 
work. 



Let us start by noticing that in numerical work is difficult to resolve the various con- 
tribution to the distribution of singular triangulations coming from the various peaks 
geometrically found by our analysis. The resolving power depends, among other param- 
eters, on the size of the triangulations, and as a rough indicator, the larger the size the 
bigger is the set of sub-dominating singular triangulations which come into play. Obvi- 
ously the first sub-dominant terms are the most relevant ones, and as suggested in the 
previous section, an interesting value to look at for comparison with Monte Carlo data is 
the value of the inverse gravitational coupling corresponding to the pseudo-critical average 
incidence (6(2, 4) s i ng ){N 4 ). As recalled there, by solving for &2 the equation 

(b(4,2) sing )(N 4 ) = 3( A ^ ) _ 1 ), (159) 

we obtain the value of /cJf^A^) corresponding to which we expect to see a clear signature 
of the dominance of singular geometries in the set of triangulated 4-spheres of volume N4. 
This is actually a pseudo-critical point, the location of which depends on JV4. Numerically 
one finds that as the volume N4 of the triangulation increases, the corresponding pseudo- 
critical point kf^iN^) increases too, (see e.g. ||). Simulations and extrapolation to 
triangulations with size N 4 = 48000 and A 4 = 64000 locate the corresponding k% rit (N A ) 
at 1.267 and 1.273, respectively. 



According to ( [L56|) the actual dependence of the number of dominating peaks, h(N4), 
as a function of the volume N4 of the triangulation, is linear according to 

h(NA = - 1, (160) 

v ; 16000 ' v ' 

for N 4 (S A ) > 32000, where the actual value, (51nft = 1/16000 and f = -1), of the 
constants comes from comparison with the numerical data provided at A^S 14 ) = 32000 



by M. With this expression of h(N 4 ) we obtain, from ( |159| ) and ( 158 ), the table ||. The 
agreement between the analytical pseudo-critical points and the Monte appears surpris- 
ingly good, and suggests that the identification of our k% with the pseudo critical 
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N 4 


h 


Analytical kf u (N 4 ) 


Monte Carlo kf tt {N 4 ) 


32000 


1 


1.25795 


1.258 


48000 


2 


1.26752 


1.267 


64000 


3 


1.27466 


1.273 



Table 3: The value of the analytical pseudo-critical points k2 nt (N 4 ) versus their Monte 
Carlo counterparts. These values are computed under the hypothesis that the linear 
dependence of h(N 4 ) from N 4 (S 4 ) is given by h(N 4 ) = iV 4 /16000 - 1. 

k% (N 4 ) found in Monte Carlo simulations is not a mere coincidence. An important im- 
plication of this identification, if correct, is that the growth with N 4 of k^ (N 4 ) is due 
to the increasing contribution of the sub-dominating singular triangulations. This result 
provides a nice explanation to the fact that Monte Carlo data seem to indicate that the 
major part of the finite size effects come from the crumpled phase 

By extrapolating the actual measurements, the Monte Carlo simulations locate the 
critical point around k\ ~ 1.327 or around k\ ~ 1.293, (depending if the data fit used 
is modeled after a second order or a first order transition, respectively)^]. Again, our 
analytical result k^" 1 — 1.3093 appears in quite a good agreement with the numerical data, 
(curiously enough our k% is, with a good approximation, the average of the above two 
numerical data), and moreover its analytical characterization provides a natural entropic 
explanation to the structure and location of the associated finite-size pseudo-critical points. 

Another distinct feature of recent numerical works concerns the bimodality in the 
distribution of singular vertices seen during Monte Carlo simulations exactly around 
k2 lt (N 4 = 32000) ~ 1.258, ||. In this connection, particularly interesting are the papers 
Q and fl3|1 , where long run histories (at N 4 (S 4 ) = 32000) provides a reliable measure- 
ment of the average maximum vertex order near the critical point. In these simulations 
the system wanders between two states characterized by two quite distinct values of the 
average maximum vertex order. In one case, this maximum is close to 3000, while for 
the other the figure is close to 1000. A correlation analysis shows that this metastability 
corresponds to tunneling back and forth from a branched polymer state (average vertex 
order ~ 1000) containing no singular vertex and a crumpled state (average vertex order 
~ 3000) with one or two singular vertices. 

According to our analysis, this behavior is the one exactly coded into the entropy 
formulae Q152| ) and ( |153j) which exactly describe a finite size tunnelling between a crumpled 
state (described by (|152|) ) and a branched polymer state (described by (|153|) . A good 
indication of the average vertex order, as we approach the transition point for increasing 
&2, is provided by ( |143|) At N 4 (S i ) = 32000 this analytic formula yields 

o 11/15 n 

(Vol(a°)) N4=32 ooo = ^5—^(7^) • ( N * = 32000) ~ 3400. (161) 
Conversely, if we approach the transition point by lowering &2 j then a reliable indication 
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is provided by (145)- Explicitly we get 



o 29/15 2 q 

(Vol(a°)) poly = -^-Exi—) ■ (N 4 = 32000) ~ 1770. (162) 

Such results appear quite in reasonable agreement with the values of (Vol(a )) at 4 =32000 
obtained during the simulations and mentioned before. Such data suggests that the bi- 
modality seen in the numerical simulation has its origin in the presence of sub-dominating 
singular triangulations. In particular, due to finite size effects the set of subdominating 
singular triangulations S 4 ,, for h = 0,1,..., 6 seems to provide a met ast able cluster of 
configurations that entropically dominate the crumpled state. 

Taken at face value, this set of results seem to indicate, at least to the indulgent reader, 
a variety of viewpoints on the actual status of a theoretical interpretation of the numerical 
simulations: 

(I) The bimodality as well as the implied first order interpretation of the transition 
between weak and strong coupling is a finite size effect related to: (i) The saturation of 
the triangulations of {§ 4 } with § 4 S in the strong coupling phase; (ii) The slow dependence 
of the specific entropy, lnc(5' 4 ; h), of {§ 4 } from the parameter h controlling the volume of 
the singular part of the triangulation. This slow variation may be responsible of the fact 
that the tunnelling does not disappear as the volume of the triangulations increases. 

Obviously, this latter remark can be easily turned inside out to favour a less optimistic 
point of view: 

(II) The slow /i-variation in In c(S 4 ; h) may well be such as to mantain the bimodality 
for larger and larger volumes: we have a genuine first order transition. 

It is rather clear that our analysis, being based on a sort of mean field approxima- 
tion, cannot distinguish clearly between such two scenarios: we need sharper entropic 
estimates. Even if shamefully low in providing answers to the headlines that numerical 
simulations score, we wish to conclude with a final example pointing to a constructive 
way of using our analytical entropy estimates. This final point concerns the &2 depen- 
dence of the two normalized cumulants of the distribution of the number of vertices of 
the triangulation, c\[N±; fo) and 02(^4; ^2) whose analytic expression is explicit provided 
by (|2"8|) and (p9|). Strictly speaking, these expressions are accurate only near the actual 
critical average incidence 60, however we can use them quite safely in a rather larger range 
of variation of &2, (due to the slow variation of 6(2,4) as a function of A^)- Accurate 
Monte Carlo measurements of such cumulants have been reported in ||, and by referring 
to these data for = 32000, the comparison between MC-data and our analytic results 
for ci(N±; fo) and C2(N±; fo) are shown in table The agreement between the analytical 
cumulant ci(k2',Ni) and its Monte Carlo counterpart is particularly good; (note that for 
a better comparison with the numerical data we have actually used in ( |29| ) an average 
between 6(4,2)1^=0 and 6(4,2)1^=1 so as to shift from fcrf — 1-3093 to a pseudo-critical 
&2 {N4) — 1.258). Slightly less impressive is the agreement between the second cumu- 
lants, but this is to be expected since near the pseudo-critical point k<j, (N4), the second 
cumulant 02(^45^2) fluctuates quite wildly. We wish to stress that such an agreement 
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r-i ( Ma • kn\ 


r^ { Mnntdrirl) 

C- I 1 VJ. \J 1 liv v. Ut b 1 


r<i( Na - kn) 


rot Mnntdnrl) 


1.240 


0.1935053 


0.18970(12) 


0.109062 


0.141(7) 


1.246 


0.1945674 


0.19150(11) 


0.1194586 


0.144(8) 


1.252 


0.1956271 


0.19399(32) 


0.1465348 


0.254(35) 


1.258 


0.1966846 


0.19712(20) 


0.3996907 


0.316(8) 


1.264 


0.1977398 


0.20052(21) 


0.1844987 


0.118(20) 


1.270 


0.1987927 


0.20085(27) 


0.1274851 


0.118(20) 



Table 4: A comparison between the analytical values and the available Monte Carlo data 
for the first two cumulants of the distribution of the number of vertices of the triangulation. 



rests both on the rigorous asymptotics (|2l|), (|23|) and on the scaling hypotheses 



and 



1, 1 



2 = ~ 



v b(k 2 



1 

V 



lim 

k 2 ->k% rit 



>b(k 2 ) b 



const., 



(163) 



(164) 



The best agreement, used in table ^j, is obtained by choosing v ~ 0.94. Eq. ( |16<: 
but a natural consequence of the vanishing of the parameter m(b) for 6(2, 4) —. 



), is nothing 
bo', whereas 

the second condition (164) rests on a less firm ground and must be considered as a working 
hypothesis to be better substantiated. 



Some of the results discussed above show that the numerical evidence pointing toward 
a first order nature of the transition can be explained in a natural geometrical framework. 
The bimodality, which has been underlined as a strong indication that the transition is 
of a first order, is well explained by the presence of entropically sub- dominating peaks 
in the distribution of singular triangulations. Similarly to what has been argued by Cat- 
terall et al. ||, the system tunnels among such distinct sub-dominant configurations 
with some of these configurations being meta-stable for A^4 finite, (especially those with 
h ~ 0, 1, . . . which dominate the crumpled phase, and those for which h >> 1 character- 
izing the branched polymer phase). Of course the analytical arguments provided by us 
are all based on a kind of mean-field approximation, since we consider only a restricted 
class of triangulations. Mean-field analysis is in general not very reliable when it comes to 
predicting the order of a phase transition. However, in this case we have seen that com- 
bined with an additional scaling assumption, we get reasonable agreement with Monte 
Carlo data for both k^N^), ci(N±) and 02(^4). This might indicate a validity beyond 
that usually provided by a mean-field approximation. 

A good test of the reliability of the geometric truncation used in the present work 
is to apply it to the more complicated system of 4d simplicial quantum gravity coupled 
to Abelian gauge fields. In that system one seemingly observe a new interesting phase 
structure IS], different from the branched polymer - crumpled phase originally reported 
in p|. 
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